arXiv:1501.00874vl [astro-ph.SR] 5 Jan 2015 


Mon. Not. R. Astron. Soc. OOP. [Tl [22l (9999) 


Printed 6 January 2015 (MN IAT^X style file v2.2) 


Mass loss on the red giant branch: the value and 
metallicity dependence of Reimers’ 77 in globular clusters 

I. McDonald 1 *, A. A. Zijlstra 1 

1 Jodrell Bank Centre for Astrophysics, Alan Turing Building, Manchester, M13 9PL, UK 


Accepted 9999 December 32. Received 9999 December 32; in original form 9999 December 32 


ABSTRACT 

The impact of metallicity on the mass-loss rate from red giant branch (RGB) stars 
is studied through its effect on the parameters of horizontal branch (HB) stars. 
The scaling factors from Reimers (1975) and Schroder & Cuntz (2005) are used to 
measure the efficiency of RGB mass loss for typical stars in 56 well-studied Galac¬ 
tic globular clusters (GCs). The median values among clusters are, respectively, 
? 7 r = 0.477 ± 0.070lQog2 and ?ysc = 0.172 ± 0.024+go23 (standard deviation and sys¬ 
tematic uncertainties, respectively). Over a factor of 200 in iron abundance, 77 varies 
by <30 per cent, thus mass-loss mechanisms on the RGB have very little metallicity 
dependence. Any remaining dependence is within the current systematic uncertain¬ 
ties on cluster ages and evolution models. The low standard deviation of 77 among 
clusters (ssl4 per cent) contrasts with the variety of HB morphologies. Since 77 incor¬ 
porates cluster age, this suggests that age accounts for the majority of the “second 
parameter problem”, and that a Reimers-like law provides a good mass-loss model. 
The remaining spread in 77 correlates with cluster mass and density, suggesting helium 
enrichment provides the third parameter explaining HB morphology of GCs. We close 
by discussing asymptotic giant branch (AGB) mass loss, finding the AGB tip lumi¬ 
nosity is better reproduced and 77 has less metallicity dependence if globular clusters 
are more co-eval than generally thought. 

Key words: stars: mass-loss — stars: horizontal branch — stars: winds, outflows — 
globular clusters: general — stars: evolution 


1 INTRODUCTION 

Mass loss is known to dominate the evolution on the asymp¬ 
totic giant branch (AGB). However, the slower evolution of 
stars with M<1 Mq means more mass is lost during the red 
giant branch (RGB) phase than on the AGB. RGB mass loss 
has been much less studied, partly because it has only mi¬ 
nor effects on RGB evolution. However, it is crucial for later 
evolution, setting the stellar temperature on the horizontal 
branch (HB) and dictating the length and path of its AGB 
evolution. 

RGB mass loss is a component of most stellar evolu¬ 
tion codes, and is generally parameterised by simple rela¬ 
tions of the stellar parameters (e.g. luminosity (L), radius 
( R ), mass ( M ), surface gravity (log g), effective temperature 
(T e ff), etc.). The two parame terisa tions most common l y used 
are those of iReimer 1 (119751 1 and ISchroder fc Cunt3 (120051 . 
hereafter SC05), which we examine in this work. However, 
use of these laws requires a fiducial reference point, parame- 
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terised by a scaling factor, 77 , which determines the constant 
of proportionality in these scaling laws. 

Calibrating 77 is important in understanding the post- 
RGB stages of low-mass stellar evolution, particularly if 77 
varies with other stellar parameters, such as stellar metal¬ 
licity. No direct variation of metallicity is included in the 
Reimers and SC05 relations, however metallicity indirectly 
affects the parameterised mass-loss rate through the conse¬ 
quent variation of stellar radius and the speed at which the 
star evolves. Mass loss from RGB stars is thought to occur 
via (magneto-)acoustic processes, hence it should be largely 
in depen dent of metallicity. If a parameterisation like those 
of|R eimersl (1 1975T) and SC05 provides is an accurate repre¬ 
sentation of the mass-loss rates of RGB stars, we can expect 
77 to be independent of metallicity too. 

RGB mass loss is notoriously difficult to measure, 
relying mainly on three tracers. Blueshifted components 
in optical and near-IR lines trace mass motions close to 
the stellar surface, but this does not necessarily translate 
into the rate at which mass i s ejected from t he st ar (e.g. 
iMcDonald fe van Loonl 120071 : iM^szaros et al.l l2009l f . Sub- 
mm CO line strengths trace material far from the star, but 
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for RG B stars is only sensitive to the handful nearest Earth 
fe.g. lGroenewegenll2014l ). Infrared excesses trace dust in the 
wind, but are hard to measure conclusively and difficult to 
translate from a dust mass-loss rate to a total mass-loss 
rate, besides which little (if any) dust production is expected 
on the RGB (iMcDonald et alll2011al ; [McDonald et al1l2012l . 
and references therein). A better way to measure mass loss 
is to obtain a difference in stellar mass between two points in 
its evolution, and use this to determine a scaling parameter, 
such as 77 . This avoids many of the systematic uncertain¬ 
ties which limit the above methods and averages out any 
short-term variations in the mass-loss rate. 

Few places offer such unique laboratories to study RGB 
mass loss as globular clusters. As comparatively simple, re¬ 
solved stellar systems, we can probe the properties of their 
individual stars with accurately known distances, metallici- 
ties, abundances and ages. The stellar mass and evolutionary 
history can easily be determined for individual stars. Mea¬ 
suring the mass of horizontal branch stars in globular clus¬ 
ters allows us to determine the 77 appropriate for ~0.8-Mq 
stars over their previous evolution. 

Determining 77 is also im portant in solving the so-calle d 
second parameter problem (lFusi Pecci fe Bellazzinil 1 1997T) . 
HB morphology of Galactic globular clusters is largely set 
by a single parameter (stellar metallicity). However, a sec¬ 
ond parameter is needed to explain the variations in mor¬ 
phological structure among clusters of a similar metallicity. 
Cluster age and helium abundanc e can both affect cluster 
morphology. iGratton et all ll2010al . hereafter GCB+10) find 
cluster age can account for most of the remaining differ¬ 
ences. If metallicity and age account for the entirety of HB 
morphology, we should expect little variation in 77 among 
different clusters, as both age and metallicity are taken into 
account in the stellar evolution models we will use. 

This manuscript is structured as follows: in Section [2] 
we collate the existing HB star masses, stellar compositions 
and cluster ages required to compute 77 . We also discuss the 
models we use to describe mass loss from RGB stars. I 11 Sec¬ 
tion [3] we describe existing and new stellar evolution models 
used to derive the mass-loss histories of each cluster’s stars. 
In Section [4] we calculate our final values of 77 and their 
uncertainties. In Section [5] we compare these results to pre¬ 
vious empirical and semi-empirical determinations, discuss 
their implications for the evolution of AGB stars in globular 
clusters, and the veracity of the variations with metallicity 
we find. 

During the course of writing this manuscript, it was 
found that publsihed stellar evolution codes were insufficient 
for our purposes. This led us to produce the new evolution¬ 
ary tracks described in Section[3] The application of all stel¬ 
lar evolution models impart poorly quantified errors to the 
results. We have quoted estimates for these errors in the 
text, but provide appendices which contain more detailed 
discussion. In Appendix m we explore the various determi¬ 
nations of cluster ages. In Appendix [B] we compare our new 
evolutionary tracks to other models. In Appendix [C] we dis¬ 
cuss various errors associated with our derivation of 77 . 


2 DATA COLLECTION AND 
HOMOGENISATION 

2.1 Data 


2.1.1 Horizontal branch masses 

GCB+10 provide the minimum, median and maximum 
masses of horizontal branch stars in 78 globular clusters, 
comprising of 74 clusters with Hubble Space Telescope im¬ 
agery, and 49 clusters for which they have ground-based ob¬ 
servations, with some overlap between the two. These 78 
clusters represent our base data set, from which we use those 
with accurately determiend ages and metallicities in pub¬ 
lished literature (59 clusters). 

Strictly, the ‘minimum’ and ‘maximum’ values of 
GCB+10 are the 5th and 95th centile masses. This equates 
to 1.65cr if one assumes a Gaussian mass distribution, though 
this is not implied by their observations. These masses 
assume that the stars have primordial helium abundance 
(GCB+10 also give the required helium enrichment to ex¬ 
plain the horizontal branch morphology at fixed stellar 
mass). They are reproduced in Table [T) 

GCB+10, in their Section 2.3, note that there is a large 
uncertainty in the maximum mass of HB stars in each clus¬ 
ter. For each star, they derive a mass based on that star’s 
colour, then apply a mass-based correction to revert that 
star to its position on the zero-age HB (ZAHB). The un¬ 
certainty in that correction, and in the amount of time the 
star has been on the HB, leads to an uncertainty in the 
ZAHB mass, which is most severe for the maximum-mass 
HB stars (their equation ( 8 )). This large uncertainty pre¬ 
cludes us from calculating accurate values for 77 for the high 
end of the mass distribution within each cluster. 

Conversely, the potential for helium enrichment among 
the stars with the minimum HB mass means that estimates 
of the evolutionary history of these stars are highly uncer¬ 
tain. This precludes us from calculating average values of 
77 for this opposite end of the HB either. For the majority 
of the following discussion, we are therefore limited to the 
median HB masses of GCB+10. 


2.1.2 Stellar compositions 


In order to obtain a reliable initial stellar mass from evo¬ 
lutionary models, we need to know the star’s composition. 
Typically, this is reduced to three parameters: the global 
metallicity, the helium abundance, and the a-element en¬ 
hancement, as these have the largest effect on the stellar 
isochrones. 

[Fe/H] is usually taken as the tracer for global metal¬ 
licity. We adopt this from the 2010 edition of the Harris 
globular cluster catal ogue (lHarrisI ll996l N Fl . Harris uses the 
metallicity scale of ICarretta et all ( 20091 ). which is based 
on more -recent, higher-resolu tio n spectroscopy than the tra - 
ditional IZinn fe West! dl984l ') or ICarretta fe Grattonl (Il997lf 


1 The unpublished documentation for this edition can be found 
at http://arxiv.org/abs/1012.3224. 
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Table 1. Data for each cluster, showing metallicities, derived ages, initial masses for the TP-AGB stars, observed current HB 
(from space- (S) and ground- (G) based observations), and mass-loss efficiencies 77 ^ and 7730 Additional systematic errors of At)r = 

and Ar]sc =lao 23 aPPty- 


Cluster 

[Fe/H] 

Age 

Minit 

Median HB mass 

77 for median stars (internal errors) 



(dex) 

(Gyr) 

(M 0 ) 

(M 0 ) 


m 




NGC 104 

-0.72 

11.95 

0.69 

0-886 

0.648 

± 

0.001 

0 452 + U - U2t * 
u.^toz _o 031 

0 1 75 + U - UU8 
6 - 1/0 _o.o 05 

NGC 1261 

1.27 

9.99 

0.50 

0.876 

0.683 

± 

0.007 

0 473 +°- 058 
-0.066 

n 1 ss + 0.020 

u .100 -0.020 

NGC 1851 

1.18 

10.07 

0.47 

0.880 1°;°22 

0.664 

± 

0.005 

0 505 + 0 052 

U.OUO -0.056 

0 192 + 0 019 

u.iy^r -0.017 

NGC 1904 

1.60 

12.14 

1.04 

0.814 ±°;“° 

0.607 

± 

0.003 

0 520 + 9,985 
u.ozu _o.o 93 

n 1 QQ +0.024 

U.18y -0.030 

NGC 2298 

1.92 

12.50 

0.77 

0.801 ±°;“g 

0.663 

± 

0.007 

0 393 +°- 074 
u.o»o _o .093 

0 143 + 0 022 
-0.029 

NGC 2419 

-2.15 

12.62 

0.93 

0-796 ±°;°g£ 

0.676 

± 

0.004 

0.4451 . 

0.158 + . 

NGC 2808 

1.14 

10.36 

0.70 

0.876 ±g;g2g 

0.608 

± 

0.000 

0 587 +°- 972 
u.ocw -0.072 

0.208 t° 0 0 0 l° 3 

NGC 288 

1.32 

11.30 

0.41 

0.843 

0.603 

± 

0.004 

0 541 +0 - 049 

U-CKfci _ 0 053 

0.193 +°;°“ 

NGC 3201 

1.59 

10.86 

0.49 

0.840 +g 

0.654 

± 

0.006 

0 492 + 0 065 
-0.071 

nisi +0.021 
U.181 -0.019 

NGC 362 

1.26 

10.05 

0.56 

0 S7 C \ +0-026 

U.OfO _0.024 

0.680 

± 

0.011 

0 475 + 0 072 
-0.077 

n 1 SS +0-024 
u.xoo _o.023 

NGC 4147 

1.80 

11.98 

0.72 

n SI 9 + 0.025 
U.812 _ 0 023 

0.648 

± 

0.006 

0 454 +°- 072 
-0.087 

0 167 + 0 021 
U- 10 ' -0.030 

NGC 4372 

-2.17 

12.73 

0.65 

0 794 + 0 026 
u. _0.025 

0.664 

± 

0.000 

0 394 +°- 047 

u.o»4: -0.072 

0 141 +0- 015 
e. 1^1 _o .028 

NGC 4590 

-2.23 

11.70 

0.64 

n 01 q +0.024 
u .010 _ Q 023 

0.704 

± 

0.001 

n qci +0.055 
u .001 - 0.086 

n 197 + 0.016 
U - 1Z ‘ -0.029 

NGC 4833 

1.85 

12.33 

0.73 

0.805 ±° 0 ° 0 ll 

0.650 

± 

0.007 

0 431 + 0 072 
_o .088 

n 1 ^s + 0.020 
6.100 -0.029 

NGC 5024 

- 2.10 

12.42 

0.62 

0.800 t° 0 ;”g 

0.657 

± 

0.001 

0 426 + 0 053 
u.^tzo -0.071 

0 153 + 0 019 
u.xoo _o.026 

NGC 5053 

-2.27 

11.88 

0.53 

0.809 t°o °o 2 2 2 2 

0.709 

± 

0.004 

0 ^95 +0-055 
u.ozo _o.o90 

0 113 +°- 015 
u .110 _o.o27 

NGC 5272 

1.50 

11.59 

0.39 

0.829 1°;““ 

0.670 

± 

0.001 

0 405 +°- 033 
U.^tUO _o .055 

n 1 C 7 +0.009 
U.IO/ _o.oi 5 

NGC 5466 

1.98 

12.53 

0.51 

0.799 

0.705 

± 

0.009 

0 9S1 +0-064 
0-Z81 -0.102 

0 104 +°- 022 

U.lU^t -0.035 

NGC 5694 

1.98 

13.64 

0.91 

0.780 h 0 0 0 ol°8 

0.644 

± 

0.006 

0 SS9 +0-073 

u.oo^ _o .092 

0 136 + 0 025 
u.xou _o.026 

NGC 5824 

1.91 

13.06 

0.52 

0 791 + 0 024 
u. < —0.024 

0.651 

± 

0.004 

0 392 + 0 045 
u.o»z, _o.o68 

0 143 +°- 012 
-0.020 

NGC 5897 

1.90 

13.00 

0.00 

0 792 + 0 015 
_ 0 .015 

0.648 

± 

0.001 

0 402 + 0 008 

U.^tuz -0.032 

0.146 t . 

NGC 5904 

1.29 

11.06 

0.62 

0 850 +0-025 
U.OOU _o 024 

0.627 

± 

0.001 

0 513 +0 054 
U .010 -0.057 

0.189 

NGC 5927 

-0.49 

11.45 

0.71 

0.910 +0 051 

0.642 

± 

0.000 

0.548 +. 

0.196 +. 

NGC 5946 

1.29 

11.63 

1.88 

0.839 +° 

0.612 

± 

0.000 

0 507 +0 - 148 

U.OU/ -0.107 

0 185 + 0 052 
6.100 -0.029 

NGC 5986 

1.59 

11.87 

0.72 

0.819 

0.611 

± 

0.001 

0 525 +0 - 061 

U.OZO -0.074 

0 190 +°- 017 

u.x^u -0.023 

NGC 6093 

1.75 

12.59 

0.74 

0.802 

0.641 

± 

0.010 

0 431 +°- 077 
u.^o± _o .092 

0 159 + 0 023 

u.xu» -0.032 

NGC 6101 

1.98 

11.85 

0.63 

0.812 

0.665 

± 

0.002 

0 433 + 0 058 

u.^too -0.076 

0 155 + 0 020 
6-100 _o .022 

NGC 6121 

1.16 

11.81 

0.50 

0.843 

0.664 

± 

0.011 

0 402 + 0,046 

U.^tuz -0.058 

0 160 +°- 013 
U.1UU -0.015 

NGC 6171 

1.02 

12.39 

0.83 

0.843 

0.650 

± 

0.001 

0 404 + 0 038 
-0.047 

0 160 + 0 010 
u.iuu -o.oii 

NGC 6205 

1.53 

12.00 

0.42 

0.820 

0.601 

± 

0.004 

n cqi +0.050 
u .001 _o .063 

0 191 + 0015 

6* 161 -0.017 

NGC 6218 

1.37 

12.51 

0.74 

0 817 +0-027 
U.81/ —0.025 

0.594 

± 

0.008 

0 501 + 0 069 

U.OU1 _o.071 

n 1 so + 0.020 

U.18Z -0.017 

NGC 6235 

1.28 

11.70 

1.88 

0 SSS +0.053 

U.8d8 _o .046 

0.617 

± 

0.051 

0 . 749 1. 

0.260 + . 

NGC 6254 

1.56 

11.78 

0.66 

O S99 +0-024 
u.ozz _0 021 

0.596 

± 

0.008 

0 557 +0-076 
u.oo/ _o .084 

n 1 QS +0-020 
6-lyo -0.024 

NGC 6266 

1.18 

11.96 

0.78 

0.839 

0.609 

± 

0.002 

0 491 + 0,055 
-0.053 

nisi + 0.017 
0.181 -0.011 

NGC 6273 

1.74 

12.28 

0.45 

0.808 

0.613 

± 

0.000 

0 512 +°- 041 

U.Oiz _o.o56 

0 185 + 0 011 
6.180 -0.018 

NGC 6284 

1.26 

11.24 

0.58 

0.848 «;S 2 24 

0.583 

± 

0.005 

0 575 +°- 071 

U.O/O - 0.068 

0 199 + 0 027 

v.i-vv -0.018 

NGC 6341 

-2.31 

12.62 

0.57 

0.795 ±0.811 

0.686 

± 

0.005 

0 346 + 0 057 

u.o<±u -0.089 

0 116 + 0016 
6.110 -0.029 

NGC 6352 

-0.64 

11.21 

0.59 

0 912 + 0 029 
-0.049 

0.667 

± 

0.007 

0.468 + 0 033 

0.180 + 0 012 

NGC 6362 

-0.99 

12.04 

0.47 

0.853 +0 ^3 

0.652 

± 

0.011 

0 419 + 0 038 
-0.046 

0 164 + 0010 
U,1U ^ -0.009 

NGC 6397 

-2.02 

12.56 

0.40 

0.798 j±°23 

0.635 

± 

0.004 

0 /LfiS +0.046 
0.468 _o.063 

0 165 + 0 016 
6.100 -0.018 

NGC 6535 

1.79 

12.11 

0.79 

0.810 j -° Q °f 4 

0.642 

± 

0.023 

0 460 + 0,120 
u.^ou - 0.132 

0 169 +°- 035 
u.ioy -0.046 

NGC 6584 

1.50 

11.37 

0.73 

0 SSS +0.025 
u .000 —0.024 

0.683 

± 

0.003 

0 390 + 0 060 

u.o»u -0.082 

0 152 + 0 017 
u.ioz -0.026 

NGC 6637 

-0.64 

11.52 

0.92 

0 905 + 0,037 
u.yuo _0 905 

0.671 

± 

0.006 

0.447 j; 0 - 047 

0.175 +°- 015 

NGC 6652 

-0.81 

11.86 

0.67 

0 S77 +0-030 

U -°‘ * -0.029 

0.650 

± 

0.001 

0 444 +0- 028 
-0.032 

n 1 7S +0-008 
U.I/O -0.006 

NGC 6681 

1.62 

12.23 

0.80 

0.812 t°o°a?2 

0.603 

± 

0.000 

fl 595 +0-063 
U.ozo -0.075 

0 190 +°- 018 

U.IUU -0.024 

NGC 6712 

1.02 

11.00 

0.00 

0.871 ±g;gi| 

0.668 

± 

0.004 

0 . 447 1. 

0.173 i. 

NGC 6723 

1.10 

12.26 

0.86 

0.839 +° ;°g° 

0.644 

± 

0.000 

0 419 + 0 041 
-0.053 

0-164 

NGC 6752 

1.54 

12.04 

0.44 

n 010 +0.018 

U .818 _o.oi 7 

0.607 

± 

0.005 

0 518 +0-052 
u .010 _o .066 

0.188 

NGC 6779 

1.98 

12.79 

0.67 

0.795 ±°;g|l 

0.652 

± 

0.001 

0.409 

0 146 + 0017 

U.l^O -0.021 

NGC 6809 

1.94 

12.49 

0.57 

0.800 

0.658 

± 

0.005 

0.408 

n 1 4S +0-017 
0.148 -0.023 

NGC 6838 

-0.78 

11.57 

0.74 

0 SS7 +0.030 
U.88/' -0.030 

0.659 

± 

0.000 

0 447 +0- 031 
-0.034 

0 174 + 0 009 
U -1'* -0.007 

NGC 6934 

1.47 

10.90 

0.59 

0 845 +°- 023 

u.o^to —0.022 

0.673 

± 

0.008 

0 440 +°- 069 
u.^u -0.078 

0 169 + 0 020 

U.xuy -0.022 

NGC 6981 

1.42 

11.02 

0.63 

0.844 +°+ 2 5 s 

0.680 

± 

0.000 

0 414 +°- 048 

-0.059 

0 162 + 0 014 
U.IOZ -0.017 

NGC 7078 

-2.37 

12.33 

0.75 

0.800 f° 0 ° 0 f 5 

0.688 

± 

0.004 

0 361 + 0 070 

u.oux -0.097 

0.115 lo’o27 

NGC 7089 

1.65 

11.59 

0.83 

0 S9S +0-026 

u.oz,o _o .023 

0.628 

± 

0.010 

0 511 +0- 094 
6.011 -0.101 

0.186 t° 0 ° 0 f 2 

NGC 7099 

-2.27 

12.83 

0.65 

0 792 + 0 025 
u.iy* _o .025 

0.665 

± 

0.001 

0 391 + 0 049 

u.o^x -0.076 

0.135 1°;°29 


masses 

+0.050 

■—0.062 
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scales (hereafter ZW84 and CG97). In practice, we can ex¬ 
pect the uncertainty of the adopted metallicities is of order 
±0.1 dex or less. We adopt this value as a global uncertainty 
on the absolute [Fe/H]. 

We make the assumption that median HB stars have 
little helium enhancement due to internal pollution, which 
implies a star-formation rate that declines relatively quickly 
from the cluster’s birth. The near-constant offset between 
the median and maximum masses in most clusters reported 
by GCB+10 suggests this is typically true. However, this 
may not be a valid assumption for a small number of clus¬ 
ters: GCB+10 single out NGC 2808 as a case in point. We 
return to this point in Section EH The uncertainty that he¬ 
lium abundance gives to the mass of HB stars is complex 
and poorly determined. We enter a detailed discussion in 
Appendix IC31 where we attribute an additional ~0.01 Mg 
uncertainty to the mass of HB stars. 

The helium fraction is usuall y taken to be a simi - 
lar formulism to that adopted by ISalaris fe Weiss! (120021 ') 
(T = 0.23 + 3 Z), which increases the global helium abun¬ 
dance only marginally from cosmological abundances, based 
on the cluster metallicity. The picture in globular clusters 
is partly clouded by the slight spread in age within each 
cluster. The latest-forming stars in a globular cluster are 
likely to be enriched with helium, nitrogen and sodium a nd 
depleted in oxygen (among other elements; iFenner et al.l 

l2004l ; lKarakas et al.ll2006l : lGratton et al.ll2010blj . These stars 

evolve faster and will not be well-represented by tradi¬ 
tional stellar isochrones, nor will helium-enriched stars have 
the same zero-age HB temperature as their primordial- 
composition counterparts, creating some uncertainty in the 
HB mass (e.g. GCB+10). Potentially large variations in he¬ 
lium abundance that may exist 24-40% ) but it is d i fficu lt 
to measure the abunda nce directly (e.g. iDupree fe Avrettl 
2013; Smith et al. 2014). For our calculations, we adopt the 
default helium enrichment of each stellar evolution model, 
neglecting internal pollution. This is discussed in more detail 
in Appendix IC3I 

The a-element (O, Si, Ca, Ti) enhancement is usually 
parameterised by the ratio |q/ Fel, w hich foll ows a known 
anti- correlati o n wit h [Fe/H] (e.gjMonacoet al. 2005). How¬ 
ever, ICarne 3 |l996l ). and later iHabgoodl (l200ll + find a con¬ 
stant value of [a/Fe] = +0.3 dex for both old and young, and 
halo a nd disc clusters. This is mirrored in iRoedieer et al.l 
(120131) . who find the same mean [a/Fe] = +0.30 dex, with a 
standard deviation of 0.11 dex. Values in their compilation 
range from 0.02 dex (NGC 6626 = M28) to 0.53 dex (NGC 
6218 = M12). Most isochrones either adopt [a/Fe] = +0.2 
or +0.4 dex (the typical range of [a/Fe]). 

The factors described in this section do not include some 
unusual abundance detai ls (e.g. fac tors of up to 10 difference 
in the Cu/Fe ratio; ISnedenl 120041 ) but include the factors 
likely to significantly alter the evolution and appearance of 
the star. 


2.1.3 Cluster ages 

Most published ages of globular clusters are on a relative 
scale, however our approach requires an absolute age. Since 
that age directly sets the initial mass, hence r/ that is derived, 
the choice of absolute age scale is an important systematic 
in the final result. 


Ultimately, the absolute age is significantly affected 
by the exact metallicity scale used, typically either that 
of CG97 or ZW84. By comparing and combining dif¬ 
ferent studies, we can provide an estimate for how ac¬ 
curate the absolute calibration of age is. We use ages 
of clusters deriv ed from scaled a verag es of the follow¬ 
ing p ub lications: Salaris fc Weissl (12002 ); de Angeli et al.l 
ll2005l) : Marin-Franch et alJ (|20 q 3 ) ; Potter et al] 1 2010 ll : 


VandenBerg et al.l (120131 ') (hereafter SW02, PA+05, 


and 

MF+09, P+10 and VBLC13). Together, they provide es¬ 
timates via 11 subtly different methods for 89 globular clus¬ 
ters. In 55 clusters, the majority of these data are derived 
from the HST ACS survey: the most-extensive homogeneous 
dataset of photometry of faint stars in globular clusters. 
Minor differences between the studies arise from different 
adopted metallicities, different methods of reddening and 
different stellar evolution codes. Petails of each study can 
be found in Appendix m 

Of the 89 clusters with ages published in the above sam¬ 
ple, 18 have the full set of 11 age measurements. They are 
typically the largest, closest, best-observed clusters, with 
the best-determined parameters: NGC 104 (47 Tuc), 362, 
1261, 3201, 4590 (M68), 5904 (M5), 6171 (M107), 6218 
(M12), 6362, 6584, 6637 (M69), 6652, 6681 (M70), 6723, 
6838 (M71), 6934, 7078 (M15) and 7099 (M30). These cover 
almost the full gamut of observed metallicities and derived 
ages. We use these clusters to calibrate the ages against each 
other. 

We scale each set of cluster ages so that the average 
age of these 18 clusters is identical in each method of each 
study. We then average the values together for each study, 
and average the studies together to create a single age. We 
arrive at a normalisation constant of an age of 11.61 Gyr for 
the average cluster. The standard deviation among clusters 
is 1.19 Gyr. We take the standard deviation among studies 
of the same cluster (typically 0.66 Gyr) as the error in age of 
each cluster, except when only one study has measured the 
age, where we take that study’s age uncertainty. Figure [T| 
shows the adopted ages for each cluster and the relative dif¬ 
ferences between the ages found by each study. The adopted 
ages are listed along with other parameters of each cluster 
in Table [TJ 

Among the different studies, the typical scatter in the 
ages is of order 4 per cent. The exceptions are for the metal- 
rich clusters, where the ages of MF+09 are noticeably (1- 
2 Gyr) older than those of SW02 and VBLC13, while the 
PA+05 and P+10 ages lie in between. VBLC13 address 
this difference in the age-metallicity gradient in detail in 
their Introduction and their section 6.1.1. They attribute it 
the sensitivity of their method to the a dopte d metallicity 
of the cluster which, being from iRutledee et alJ (Il997l ). in¬ 
volves some extrapolation at the metal-rich end. We have 
maintained a straight average of the five studies for our re¬ 
sults, accepting the increased uncertainty for the clusters 
with metallicities of [Fe/H] > -1.1 dex. 

A compilation of published ages for 38 clusters, 
IRoedieer et al.1 (1201314 . was published during the prepara¬ 
tion of this paper. They give ages which are relatively similar 
to ours: they are on average 0.66 Gyr older than ours, and 
the differences have a standard deviation of 0.88 Gyr. The 
most devian t ages are those with t he highest uncertainties, 
as noted by [Roediger et al.l ll2013l) . Excluding the outliers 
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[Fe/H] 



[Fe/H] 

Figure 1. Adopted ages of globular clusters on an absolute scale (top panel) with the relative deviations in age of each included study 
(bottom panel). The large, black squares with error bars denote values and uncertainties adopted; smaller, coloured points denote the 
individual studies: SW02 - magenta asterisks; DA+05 - blue crosses; MF+09 - green plus signs; D+10 - cyan open squares; and VBLC13 
- red dots. 


NGC 5946, 6235, 6342 and 6544, their values are 0.95 Gyr 
older with a standard deviation of 0.50 Gyr. We retain this 
set of ages as a separate comparison during our analysis. 


2.2 Mass-loss models 

While a variety of formulisms have been employed to gener- 
ate mass loss f rom stars, by far the most enduring is that of 
iReimersI (|l975i ): 

T R 

M = 4 x 10~ 13 ?7r — [M© yr _1 ], (1) 

where M is the mass-loss rate, in solar units, and L is the 
luminosity, R the radius and M the mass of the star in solar 
units. The factor r/ is a scaling parameter that describes the 
efficiency with which mass is lost. 

Many attempts have been made to improve on Reimers’ 
formulism. Most have some restrictions which make them in¬ 
applicable, either: (a) they do not cover the low-mass regime 
we are interested in; (b) they are tuned to main sequence 
stars, and do not reproduce the dominant RGB mass loss 
well; (c) they do not separate AGB stars undergoing dust- 
driven winds from dustless RGB stars; or (d) they focus 
solely upon dust-driven AG B winds (examples can b e found 
in the comparative work bv ISchroder fc Cuntal2007l ). These 


later stages are not relevant to deriving r/B,. However, SC05 
covers the regime of globular cluster RGB stars, and is hence 
relevant to our current study. 

SC05 take what they describe as a “physical approach”, 
adapting Reimers’ model based on scaling laws of chromo- 
spherically driven winds. They derive the following formula: 



where g is the stellar surface gravity. In SC05, these authors 
also derive r /sc for the two globular clusters NGC 5904 and 
5927 (see Section fSTTl) . fitting r ]sc ~ 0.2, qualitatively sug¬ 
gesting that th is provides a better fit t han Reimers’ law. 
Further work in lSchroder fe Cuntd d2007l) suggests that this 
law also provides a better fit of observed mass-loss rates over 
a wide variety of masses and metallicities. By extending this 
to a wider sample of globular clusters, we can determine 
whether this remains true over considerably lower masses 
and lower metallicities. 

Though there is no physical reason why r] cannot be 
time-variant, the laws are defined such that r/ should be 
constant over time. Making this assumption, one can calcu- 
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late r] from AM by solving the time integral of the above Table 2. Adopted elemental abundances for the mesa models, 
equations, namely: All other elements set to [X/Fe] = 0. 


AMr = 4 x 10 13 77 r 


and: 


r* l 15 


r 0 M(T/T e K) 


rd t [M©]. 


(3) 


AMsc = 4 x 10 rise 


( TeS f 

V4000K/ 


1 + 


o M(T ea /T Q ) 2 

L 

4300 M(T e f[/T®)‘ 


dt [M©]. (4) 


Here, R has been substituted for VL/T p 2 ff] where T e « 
is the stellar effective temperature. The difference in mass, 
AM, can be numerically integrated for a specific g using 
stellar evolutionary tracks, allowing comparison to the ob¬ 
served change in mass, AM n b 8 , derived from the above ob¬ 
servations. 


Abundance ratio 

Adopted value 

Y 

0.2485 + 1.78 Z 

[Z/ Fe] 

+0.226 dex 

[C/Fe] 

-0.4 dex 

[N/Fe] 

+0.9 dex 

[O/Fe] 

+0.2 dex 

[Ne/Fe] 

+0.3 dex 

[Na/Fe] 

+0.3 dex 

[Mg/Fe] 

+0.4 dex 

[Al/Fe] 

+0.4 dex 

[Si/Fe] 

+0.4 dex 

[S/Fe] 

+0.3 dex 

[Ar/Fe] 

+0.3 dex 

[Ca/Fe] 

+0.3 dex 

[Ti/Fe] 

+0.2 dex 

[Cr/Fe] 

+0.1 dex 

[Ni/Fe] 

0 dex 


3 STELLAR EVOLUTION MODELS 

With the cluster parameters known, we now need to use 
stellar evolutionary models to calculate the total mass loss 
experienced to date by stars currently on the HB. 


3.1 Published models 

Stellar evolutionary models can provide initial masses, but 
each comes with its own restrictions. It is important for our 
analysis that the evolutionary tracks correctly reproduce the 
stellar parameters (L, R, T e g, g and M) and their time evo¬ 
lution. Most evolutionary tracks incorporate varying initial 
stellar mass and metallicity. However, not all allow enhance¬ 
ment of a elements which, in cool stars, primarily affects the 
atmospheric opacity, hence R , g and T e g. A proper prescrip¬ 
tion of stellar mass loss on the RGB is also required, as this 
affects gravity, causing repercussive effects in T e ff and R. 
The internal pressure is also altered such that L changes 
and, along with it, the timing of the helium flash that ter¬ 
minates the RGB. 

No publicly available set of stellar models incorporates 
both a sufficiently variable [a/Fe] ratio and a proper RGB 
mass-loss treatment. Of those t hat come close, the 2012 ver¬ 
sion of the Dartmouth models jDotter et al.ll2008l f incorpo¬ 
rates a variable [a/Fe], but does not apply mass loss while 
the star is on the RGB (A. Potter, priv. com m.). Version 1.1 
of the PARSEC models dBressan et al.l I2012T) includes mass 
loss in the form of a variable Reimers’ ig, but is restricted 
to stars of solar-scaled composition and publicly presents 
isochrones rather than evolutionary tracks. We compare to 
these and other publicly available model sets, the details 
of which are given in Appendix [B] However, for a proper 
derivation of r;, we must create our own set of models which 
incorporate both these effects. 


3.2 Mesa evolutionary tracks 


We have created a grid of evolutionary tracks using 
the MESA (Modules for Experiments in Stellar Astrophysics) 


stellar evoulution code l|Paxt,on et al .Ifioul . 1201 .'ll ). This grid 
spans [Fe/H] = -2.6 to -0.2 dex in 0.2 dex intervals, and ini¬ 
tial stellar masses of 0.76 to 0.94 M© at 0.02 M© intervals. 
The grid is not complete, but models are run to cover the 
entire range of globular cluster star initial masses and metal- 
licities, including the range of likely uncertainties. Models 
were run for 15 Gyr of evolution, but generally only main¬ 
tain convergence as far as the last or second-last thermal 
pulse. All models are complete at least to the start of the 
thermally pulsating AGB. 

Default MESA parameters were assumed in most cases, 
with the following exceptions. Helium abundance was set at 
Y = 0.2485+1.78Z and Z© was taken as 0.0152. Metals were 
fractioned to impart [q/Fel ~ +0.3 dex, but with detailed 
abundances based on iRoediger et ahl (|2013l l and references 
therein, as listed in Table [5] The initial pre-main-sequence 
models began with an initial core temperature of 5 x 10 s 
K and were relaxed for 200 steps before further processing. 
‘Type 2’ opacities were used, and the kappa_f ile_pref ix op¬ 
tion was set to gs98_aF_p3. re presenting an [a/Fe] = +0.3 
dex dGrevesse fe Sau vall 199 8h. Henye y mixing-length the¬ 
ory was adopted ( Henvev et all 119651 1 with omlt = 1.92 
to match the RGB of 47 Tucanae (see Appendix IB2I) . The 
which_atm_option atmospheric boundary condition was set 
to photosphere_tables. 

Mass loss can be included into MESA in a variety of 
built-in parameterisations. We applied Reimers’ mass-loss 
law to the entire stellar evolution. We compute two sets 
of models, for ??r = 0.4 and ??r = 0.5. The law of SC05 
is not a presently available mass-loss prescription in MESA. 
The inherent differences in stellar evolution between a star 
undergoing Reimers-like mass loss and a SC05-like mass loss 
impart an additional uncertainty to our determination of 
Vsc- 


A detailed discussion of the systematic uncertainties im¬ 
parted by the choice of stellar evolution model, stellar com¬ 
position and stellar physics can be found in Appendix [C] A 
comparison of the MESA tracks to the Dartmouth and PAR¬ 
SEC models can be found in Appendix IB2I 
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A calibration of RGB mass-loss efficiency 7 



Figure 2. Lines of constant mass in metallicity—age space from 
the MESA models. The literature-derived locations of globular clus¬ 
ters are shown as black points. 
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Figure 3. The error in initial mass arising from choice of stel¬ 
lar model, including contributions from convective overshooting, 
[ft /Fe] and natal helium abundance, as described in Appendix [Cl 



3.3 Initial stellar masses 

Figure [5] shows the initial masses of RGB tip stars as taken 
from our MESA models and the Dartmouth isochrones at 
[cr/Fe] = +0.3 dex (linearly interpolated between their +0.2 
and +0.4 dex results). Overlain are the globular clusters. 
Most clusters have RGB tip stars with initial masses be¬ 
tween 0.76 and 0.90 M@. The most-metal-rich have higher 
masses, but still below 1 Mq . Only the three youngest clus¬ 
ters (in order of increasing age, Palomar 1, Terzan 7 and 
Palomar 12) have masses that approach or exceed 1 Mq. 
We do not model these younger clusters in this work. 

Uncertainties in the derivation of initial mass arise from 
uncertainties in the clusters’ metallicities (A[Fe/H] = 0.1 
dex; Section [2.1.211 . ages (Section 12.1.311 and stellar proper¬ 
ties. For the stellar properties, we consider uncertainties due 
to the choice of stellar evolution code, helium and a-element 
abundances, and convective overshooting as described in Ap¬ 
pendix [C] Figure [3] shows the associated error arising from 
the parameters. Typically this combined error is ~0.01-0.03 
Mq, rising with metallicity. 


4 MASS-LOSS EFFICIENCY ON THE RGB 
4.1 Deriving 77 

In this Section, we fit 77 r and 17 sc to match the mass lost by 
HB stars in each cluster between their birth (Section 13.311 
and their arrival on the ZAHB (GCB+10; Section 12.1.11) . 
This calculation is performed as follows. For each cluster, we 
identify the four ? 7 r = 0.40 MESA models which bracket it in 
ZAHB age and [Fe/H], For each model, we identify the mass 
lost between the star’s birth and the ZAHB. Interpolating 
across the four models in age and [Fe/H], we find the mass 
loss appropriate for that cluster’s age and [Fe/H]. We repeat 
this for the 77 R = 0.50 models. We construct a linear fit 
between the mass lost in the 77 r = 0.40 and 0.50 models, and 
identify where the observed mass loss falls on that linear fit, 
giving us an 77 for that cluster. 

This process is repeated for points at the extremities of 
the age and metallicity error range (Section 13.31) to provide 
an error in ? 7 r which is relative compared to the values for 
the other clusters in this study. An additional systematic er¬ 
ror applies, arising from the choice of stellar evolution model 
and its parameters (also Section [+3]), which is quoted sepa¬ 
rately. 

The analysis described in this section is performed again 
using the law of SC05. The ? 7 sc for each stellar evolution 
model is calculated as the value that provides the same 
amount of mass loss from the star’s birth to the ZAHB. The 
calculation is then performed in the same way, although an 
additional error is added to account for the uncertainty cre¬ 
ated by adopting r/sc rather than ? 7 r (Appendix IC4I) . 

Figure [4] shows our derived values of 77 for the median 
stars in each cluster. 


4.2 Application to the mimimum and maximum 
helium abundance 

As noted in Section [ 2 . 1.11 GCB+10 provide a range of HB 
star masses for each cluster. However, the applicability of 
this range here is limited due to the uncertain helium abun¬ 
dance. 

GCB+10 note a small (0.03-0.04 Mq), near-constant 
offset between the median and maximum masses in almost 
all cases, which translates to a change in 77 R of 0 . 1 , and a 
change in rjsc of 0.03. This would be indicative of a spread in 
77 within each cluster of <t( 77 r) = ±0.05 and <t(t 7 r) = ±0.015. 
However, the correction applied by GCB+10 is too uncertain 
to say whether this represents an intrinsic spread in 77 which 
might arise from a lower helium enrichment. 

Conversely, the maximum 77 can be a useful measure 
if the helium enrichment in that cluster is low. Anomalous 
enrichments of AY<0.01 are seen by many authors (e.g. 
Catalan et al.l l2009inVillanova et al.l 120091 : IValcarce et al.l 
2O13e0) but severe enhancements of up to Y ~ 0.4 are 
estimated by other authors in the most massive clusters 
(e.g. NGC 2808, NGC 6441, c o Cem ld’Antona et al.1120051 : 
iPiotto et al.ll2005l : ICaloi fe D’Antonall2007l ~). 

Comparison of mesa tracks of 0.90 Mg stars at [Fe/H] 
= -0.60 dex show that a helium enrichment of AY = 0.005 


2 However, note hig h er values of helium enr ichment in M4: 
IVillanova et al.1 ll2012h ; iDalessandro et al.l d2013l '> . 
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[Fe/H] 


Figure 4. Top panel: initial masses of HB stars in each globular cluster (upper, red points with solid error bars) with the current median 
masses of those stars, from GCB+10 (lower, blue points, having negligible error bars). Green lines give the canonical white dwarf mass of 
0.53 d= 0.02 Mq. Central and bottom panel: variations of the associated median mass-loss efficiencies ( 77 ^ and t/sc)- Red points show the 
observed values per cluster and their associated uncertainties under the assumption that there is no helium enrichment. Green, dashed 
lines show the average value (thick line) and standard deviation (thin lines). Blue, dotted lines show a linear fit (thick line) along with 
the standard deviation from that fit (parallel thin lines) and error in the slope (other thin lines). Errors on clusters with [Fe/H] > -0.64 
are not always fully defined. 
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A calibration of RGB mass-loss efficiency 9 


will decrease in the timing of the RGB tip by 0.54 Gyr. 
This would decrease the initial mass of stars which today 
populate the HB of 0.011 Mg (full details are presented in 
Appendix IC2. 111 . This translates into a decrease of A 77 R ss 
0.022 at [Fe/H] = -1 dex and Ar/R ~ 0.033 at [Fe/H] = -2 
dex. 

GCB+10 model a large scatter in the difference between 
the median and minimum masses of HB stars, but the aver¬ 
age across all metallicities is ~0.057 M@. If the change in HB 
mass is caused entirely be helium enrichment, this translates 
to an increase of AT ~ 0.026 between a cluster’s median- 
mass and minimum-mass HB stars. Since GCB+10 model 
that 68 per cent of clusters have a maximum Y of 0.293 or 
less, this is not an unreasonable result. Alternatively, if the 
change in HB mass is caused entirely by RGB mass-loss effi¬ 
ciency, this translates to an increase of 77 R = 0.11 at [Fe/H] 
= -1 dex and 77 r = 0.16 at [Fe/H] = -2 dex. This increase 
would be considerably greater than the spread between clus¬ 
ters. 

Without accurate measures of helium abundance in 
globular clusters, it is not possible to differentiate between 
these two scenarios, hence we only further consider the 77 
that can be applied to each cluster’s median-mass HB stars. 


4.3 Final average values of 77 

From the data presented in Figure [T] and Table |T] we derive 
an average value across all clusters of r/n = 0.477 ± 0.070 
and ? 7 sc = 0.172 ± 0.024, where the quoted uncertainties 
refer to the standard deviation of best-estimate values for 
each cluster. However, the error budget for each cluster 
is dominated by the systematic uncertainties in determin¬ 
ing the initial stellar mass and stellar evolution pathways, 
which derive from the uncertainties in the cluster ages, he¬ 
lium content and choice of isochrones. We adopt a system¬ 
atic uncertainty in age of ±0.7 Gyr, conservatively based 
on the standard deviation of ages derived for each cluster. 
Although we use similar sources, thi s can be com pa red to 
the difference in average age between iRoediger et al.l (l2014l 'l 
and this work of 0.28 Gyr. This imparts a systematic un¬ 
certainty of At?r = 0.033 and Ar/sc = 0 .012. Our estimated 
correction for helium enrichment adds a further downward 
uncertainty of A77R = 0.037 and A?;sc = 0.013. The typical 
systematic uncertainty derived from the choice of isochrones 
is A?7r = 0.037 and Agsc = 0.014. A final, additional er¬ 
ror of Ar/sc = 0.003 is applied to account for the fact we 
do not calculate our RGB isochrones using this mass-loss 
prescription (Appendix IC4I) . Combining the systematic un¬ 
certainties in quadrature leads us to our final median values 
of 77 , which are: 

t?r = 0.477 ± 0.070tg;$g 
Vsc = 0.172 ± 0.024±]/q23, 

where we quote the standard deviation among clusters and 
the global systematic uncertainty, respectively. 

There is little noticeable variation of 17 among the Milky 
Way globular clusters, despite covering a significant range 
in parameter space, particularly in metallicity. However, 
we note a significantly larger spread in 17 in the metal- 
intermediate clusters and a tail off to lower 77 in the metal- 
poor clusters. We discuss these variations in Section Em 


5 DISCUSSION 

5.1 Comparison with previous derivations of 77 


Our derived values of ? 7 r = 0.477 ± 0.070lg;Qg2 and 
? 7 sc = 0.172 ± 0.0241 q q 23 are in good agreement with pre¬ 
vious literature. SC05 derive both 77 R and r/sc for the clus¬ 
ters M5 (NGC 5904) and NGC 5927. Their results are sum¬ 
marised in Table [3] Both clusters’ values for 77 R and 773 c are 
consistent within the combined err or budgets. Further com - 
parison can be made to field stars: ICranmer fe Saail (1201 ll i 
model 47 Galactic stars, finding 773 c = 0.2125 (with a pre¬ 
sumed uncertainty of 0.0125): also consistent within the 
combined error budget. 

The precise value of both ? 7 r and r/sc one derives de¬ 
pends on the treatment of the uppermost parts of the red gi¬ 
ant branch, where mass-loss rates reach ~2 xlO - ' Mg yr _1 . 
We have tried to incorporate reasonable errors by compar¬ 
ing several stellar evolution codes in our analysis, however 
it is not clear which stellar evolution model SC05 use. 

Discussions brought up during the refereeing process 
of this work have highlighted the need for a proper treat¬ 
ment of stellar mass loss within the comparison code. While 
the uncertainties due to variation of giant branch mass loss 
within the stellar evolution codes are incorporated into our 
systematic error, their precise treatment can alter the value 
of 77 significantly. For comparison, we have repeated our 
derivation of 17 with the Dartmouth isochrones, which do 
not include mass loss while the stars are on the RGB. The 
Dartmouth-based 77 is calculated as the substantially higher 
? 7 r = 0.550 ±0.062 and 773 c = 0.261 ±0.024 due to its hotter 
RGB tip and altered RGB lifetime (see Appendix IB2I) . 


5.2 Empirical calibration of mass loss 


The mass of a star can be calculated from observables using 
the formula: 


M = 


gL 

AnaTGG 


(5) 


For measurements of individual stars, fractional uncertain¬ 
ties in both L and T 4 can be as little as ~3 per cent. How¬ 
ever, purely spectroscopic derivations of g are normally un¬ 
certain by a factor of three or more. This make it difficult to 
achieve sufficiently accurate masses for individual stars, par¬ 
ticularly absolutely calibrated masses. However, computing 
the difference in mass between two populations of stars from 
the same observation is possible. This method has the ad¬ 
vantage of being largely insensitive to the helium abundance 
of the star. 

In collaboration with C. I. Johnson, we have previously 
adopted this approach to measure the difference between 
early AGB stars and RGB stars of similar luminosity in w 
Centauri, measuring a 26 ± 4 per cent decre ase in mass 
from the RGB to AGB (iMcDonald et al1l201ld l. Assuming 
a 12-Gyr-old cluster at [Fe/H] = -1.62, this equates to 0.21 
± 0.03 Mg, which compares favourably with the 0.23 Mg 
expected using our fit to Reimers’ law with r/n = 0.477. 

Similarly favourable results can be found in 47 Tuc, by 
co mbining the photometric temperatures and luminosities 
of [McDonald et alj (l2011bl ) with the spectroscopic log(g) 
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Table 3. Comparative values of g from SC05. Their errors are implied from their text. Systematic errors have not been applied to our 
data here. 


Cluster 

SC05 

This work 


Vn 

3SC 

VR 

3SC 

NGC 5904 

NGC 5927 

O O 

-H -H 

CO lO 

O O 

0.2 (~ ±0.04) 
0.2 (~ ±0.04) 

0-545+^“ 

0-548t°;°” 

0 194 +uu ^ 

- 0 . 016 

0 196 -1 " 0018 



[Fe/H] 


Figure 5. As Figure [4] Here, the addition of hollow magenta 
circles show the median masses predicted for stars reaching 1000 
Lq on the AGB, where dust production is expected to begin. 
Errors have been omitted for clarity, but are typically 0.025 Mg 
on the initial masses, and 0.01 Mg on the HB and at 1000 Lq. 


measurements of RGB stars by ICordero et al.l ( 2014 ) and 
of AGB stars by Johnson et al. (submitted). Although the 
errors from these more-direct measurements are still consid¬ 
erable, they constrain the possible range of ijb. to between 
~0.30 to ~0.60, re-enforcing the values we derive here. 


5.3 Implied evolution along the AGB 

The r] we have calculated above should also apply in later 
evolution, as the stars are thought to continue to lose mass 
via the same mechanisms well into their AGB evolution. In 
this section, we follow our stellar evolution models onto the 
AGB and compare their predictions to observations of AGB 
stars in globular clusters. 


5.3.1 Applying Reimers’ mass-loss law to the AGB 

Highly evolved stars are thought to lose mass primarily 
through pulsation-enhanced, dust-driven winds. Pulsations 
in the stellar atmosphere levitate material which can con¬ 
dense into dust; radiation pressure on this dust forces it from 
the star. However, the criteria for this to be the primary ejec¬ 
tion mechanism are poorly determined, partly depending on 
the strength of the pulsation and the opacity of the dust (e.g. 
IWoitkell2006l : iBladh et fill 1 20131) . In the early stages of this 
regime, where pulsation and dust production are visible but 
yet to become effective, we may still expect a Reimers-like 
law to be effective. However, if pulsation and dust driving 
enhance mass loss, the AGB should be truncated at a lower 
luminosity. 


This process should be metallicity dependent. Metal- 
poor stars are hotter, typically exhibiting weaker, shorter- 
period, less effective pulsations. They also suffer from very 
low dust-to-gas ratios, meaning dust driving is less effective 
at accelerating all the ejected material. We can therefore 
expect that the AGBs of metal-rich clusters would be more 
curtailed than metal-poor clusters. 

Dust production in globular clusters typically b egins 
somewhere between 700 Lq and 1500 Lq dMcDonald et al.l 
2009|; BweretaL 20091 ; [McDonald et al.lT2011bl lalldl ; see also 


Momanv et al.l 20121 ). It is becoming clear that the average 


star starts to develop a dusty wind close to 1000 Lq. Our 
mass-loss efficiency can be used to predict the mass of a star 
reaching this 1000 Lq point. By comparing this to the final 
white dwarf mass, we can determine how much more mass 
it should lose in a dusty wind. 

However, accurate w hite dwarf masses fo r globular clus¬ 
ters are difficult to find. I Richer et alj (Il997l ) measured the 
mass of white dwarfs in M4 (NGC 6121) t o be 0.51 ± 0.03 
Mq, a mass later refined bv IKalirai et al.l f2009) to 0.53 ± 
0.01 Mq. iMoehler et al.l d2004l) measured masses for white 
dwarfs in NGC 6752 at 0.53 ± 0.02 Mq. The initial-final 
mass function for stars at these masses is expected to be 
relatively flat (IKalirai et al.ll2009l ; iGesicki et al.ll2014l) . 


5.3.2 Results from stellar evolution modelling 

Finding the stellar mass at 1000 Lq becomes a trivial matter 
of following each stellar evolution model until it reaches 1000 
Lq, interpolating across both the rpi = 0.4 and 0.5 tracks 
to find a mass appropriate to any particular value of 77 . To 
find the mass lost in a dusty wind, we also calculate the 
envelope mass of an AGB star at the luminosity of the RGB 
tip. These are presented in Figure [5] 

We also estimate the final AGB luminosity, avoiding 
very short extensions to high luminosity during AGB ther¬ 
mal pulses. However, the lack of convergence during the fi¬ 
nal or second-to-final thermal pulse means that we can only 
give these approximately. Envelope masses at the final model 
step are typically ^0.01 Mq, which will be lost in approxi¬ 
mately ~30 000 years. Our computation of the time spent on 
the dust-producing AGB (above 1000 Lq; Figure^ and the 
time spent as an AGB star with a luminosity greater than 
that attained on the RGB, will therefore be underestimated 
by ~30 000 years. Stars exceeding the RGB-tip luminosity 
are easily identifiable as AGB stars, as they stand brighter 
than the RGB tip in infrared colour-magnitude diagrams. 

Clusters with metallicities above [Fe/H] = -1 dex are 
not shown in Figure [5] as interpolation becomes difficult at 
higher metallicities. We stress these results are for the stars 
which have had median ZAHB mass; more- or less-massive 
stars will experience more or less dusty mass loss. Results 
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Figure 6. The predicted lifetime of stars (colour scale) on the dust-producing AGB (L > 1000 Lq) and as AGB stars visible above the 
RGB tip, for the two modelled values of r /. Globular cluster initial masses, derived here, are shown as small black points. Grey colours 
indicate the star does not evolve as far as the indicated phase, clusters in yellow boxes may have some observable dusty/luminous AGB 
stars, clusters in darker (purple, blue or green boxes) should have observable dusty/luminous stars. Models were not computed in the 
hatched area, as evolution to the RGB tip takes longer than a Hubble time. 


from the mass-loss law of SC05 are not shown, as they have 
no clear terminus. 


ing nearly 5000 L n (IMcDonald et al.ll20lTbl ; lLebzelter et all 


120141 : [McDonald fe Ziilstral 2014 ). 


5.3.3 Comparison to observed AGB tip luminosities 

Figures [5] and [ 6 ] imply that the median stars in clusters have 
between 0.01 and 0.15 M 0 to lose as dusty winds, with a 
typical value being closer to 0.10 M@. Metal-poor clusters 
would produce proportionally more mass in dusty winds, 
although the presumably lower dust-to-gas ratios of those 
winds would mean less dust overall. Consequently, stars in 
metal-poor clusters should reach a higher luminosity on the 
AGB, producing dust for around ~1.5 Myr and remaining 
above the RGB tip for ~1 Myr. Conversely stars in metal- 
rich clusters should produce dust for <1 Myr. I 11 many 
metal-rich clusters, there should not be a star above the 
AGB tip for much of the time. 

Observationally, however, M15 (the most-metal-poor 
cluster) has an A GB which terminates at a luminosity be - 
low the RGB tip dBover et al.ll2006l : iMcDonald et ~al]|2010h . 
It is likely that a lack of bright AGB stars is typical 
of most or all metal-poor clusters. AGB stars in clus¬ 
ters appear to just reach the RG B tip in ai Cen ([Fe/H] 
~ - 1.6 dex; IMcDonald et al l 120091 ) and noticably exceed 
the RGB tip in luminosity by NGC 362 ([Fe/H] = -1.16 
dex: iBover et al.l 120091 ). In 47 Tuc ([Fe/H] = -0.72 dex), 
stars spend ~250 000 years above the RGB tip, reach- 


We therefore have a discrepancy whereby the prediction 
using Reimers’ law is that metal-rich clusters have no bright, 
dusty AGB stars and metal-poor clusters have bright, dusty 
AGB stars, whereas the opposite is observed. The problem 
is compounded when one observes that metal-poor clusters 
may have a slightly lower 77 R, than metal-rich clusters, lead¬ 
ing to longer-lived stars and a higher-luminosity AGB ter¬ 
minus (see Figure [4]). A departure from a Reimers’-like law 
to a pulsation-enhanced, dust-driven wind would exacerbate 
this, therefore we do not think it likely that pulsation and 
dust production greatly affect the mass-loss rates of even 
the brightest globular cluster stars. 

Adopting the SC05 law partially alleviates this discrep¬ 
ancy. The lack of a full treatment of the SC05 law means 
that the results are marginally less secure, however we can 
approximately identify the AGB terminus under the SC05 
law by identifying when a star undergoing mass loss fol¬ 
lowing the SC05 law will have lost its entire envelope. For 
metal-poor clusters ([Fe/H] < -2 dex) the AGB terminates 
around 900 000 years earlier for the SC05 law; for metal-rich 
clusters ([Fe/H] > -1 dex) it is about 400 000 years earlier. 
This is not enough to reverse the trend in AGB luminosity 
seen in Figure [ 6 ] 
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5.4 The underlying cause of variations in 77 among 
clusters 

The variance of 77 among clusters is surprisingly small: the 
standard deviation among clusters is only 14 per cent. This 
is despite their variation by over 20 per cent in initial mass 
and two orders of magnitude in metallicity, with RGB-tip 
stars that are up to 1300 K warmer and a third fainter in 
metal-poor clusters than metal-rich clusters. That this works 
so well is a testament to the accuracy of both mass-loss 
formulisms studied here. 

SC05 claim their formula should provide better accu¬ 
racy than Reimers’ due to the incorporation of physical 
princples, rather than a simple empirical fit. We note that 
both values for 77 come with a 14 per cent standard deviation 
in this work, thus one would not appear any more accurate 
than the other. However, the median internal (cluster-to- 
cluster) error we derive for 77 is also 14 per cent. This in¬ 
ternal error is driven largely by the uncertainty in cluster 
age for metal-poor clusters, and a combination of age and 
metallicity for metal-rich clusters. It is not clear from our 
data that the SC05 law necessarily provides a much better 
fit than Reimers’. 

Despite the relatively closeness in 77 among clusters, 
there are still some trends that can be seen in Figure U 
Firstly, there is a much greater variation in 77 for the metal- 
intermediate clusters (-1.7 < [Fe/H] < -1.1 dex). Outside 
of this range, there appear noticable gradients in the results 
with 77 declining at the lowest metallicities. These are re¬ 
flected in the predicted lifetimes of AGB stars in the dusty 
and super-RGB-tip phases described in Section 15.31 where 
metal-poor clusters have AGB tips higher than typically ob¬ 
served and metal-intermediate clusters have a range of ex¬ 
pected AGB termini. 


5-4-1 The gradient in 77 

To reconcile our Reimers’ 77 calculations with observed AGB 
populations for the metal-poor clusters, one or more adap¬ 
tations could be made: 


• A mass-loss law with a stronger temperature depen¬ 
dence could be used, possibly invoking a mechanism which 
more strongly affec ts hotter HB or early-AGB stars (e.g. 
iDupree et al. 120091 1. The temperature dependence would 
have to be stronger than that of SC05. 

• Slower evolution in metal-poor AGB stars in the ~500- 
2000 Lq range than predicted, allowing the wind to be lost 
over a greater amount of time. 

• The masses of the horizontal branch stars in metal-poor 
clusters could have been over-estimated by GCB+10. Figure 
[4] shows a noticable upturn in the HB mass near [Fe/H] = 

1.8 dex. If these stars can be made less massive, not only 
do the 77 of these stars become similar to that of the other 
clusters: the AGB tip will also decrease in luminosity to the 
observed value. 

• The ages of the oldest globular clusters could be made 
younger. This would increase the initial mass of stars, in¬ 
creasing 77 to values more consistent with metal-rich clusters 
and allowing more AGB mass loss to take place, truncating 
the AGB at the observed values. 


The first two possibilities are relatively unlikely. HB mass 
loss would have to increase above predictions by a factor of 
~4 to prevent super-RGB-tip stars to become visible. Slower 
evolution would change the luminosity function around the 
RGB tip in younger environments (e.g. dwarf galaxies), but 
there is no noted drop in the number stars per magnitude 
over the AGB tip in these environments. 

GCB+10 quote errors of only a few x 0.001 M@ on their 
HB star masses. In practice, however, this neglects contri¬ 
butions from the uncertain metallicity of the clusters, and 
the stellar evolution codes they use to determine H B mass. 
The difference in metallicity quoted by GCB+10 and [Harris 
(|20ld l for metal-poor clusters is ~0.03 dex. A 0.1-dex un¬ 
certainty in an individual cluster’s metallicity changes the 
HB mass by ~0.01 Mg, and ( 77 R | r/sc) by ~(0.03 | 0.01). 
While this would remove a significant part of the gradient 
in the 77 versus [Fe/H] relation, without a systematic change 
in many clusters, global [Fe/H] is not likely to be the under¬ 
lying cause of the differences. Similarly, the contribution by 
helium is not expected to be significant (Appendix IC3I) . 

The absolute choice of metallicity (Z), however, may 
be an is sue. GC B+10 use t he metal-poor Pisa evolutionary 
models llCariulo et al .lEffiOlli . which adopt \Z/X\® — 0.023, 
compared to the modern value of Z = 0.0152 ± ~0.0006 
(see Appendix IC2.2I) . This amounts to a «0.083 dex dif¬ 
ference in the zero (solar) point for metallicity, depending 
slightly on the helium abundance adopted. It is not imme¬ 
diately obvious how correction was made by GCB+10 from 
the s olar-scale d Pis a mo dels to the a-enhanced globular clus¬ 
ters. jCariulo]^^] (120041 4 provide a method for doing so, but 
it involves the adoption of a particular [a/Fe], which is not 
quoted by GCB+10. Given the likely uncertainty in the av¬ 
erage [a/Fe] adopted (~0.1 dex), one c an presume a f urth er 
~0.06 dex uncertainty in Z , based on ICariulo et al.l ([20041 . 
their section 2). The combination of these uncertainties is 
sufficient to alter the derived horizontal branch masses sig¬ 
nificantly (we estimate by up to ~0.02 M@ for the metal-rich 
clusters). However, this alters the metal-rich end more than 
the metal-poor end, and increases rather than decreases the 
gradient with metallicity for both 77 R and r\sc versus [Fe/H], 

The choice of evolution model is also significant. The 
Pisa HB models for Z = 0.0002, M = 0.65 Mq return a 
ZAHB temperature of ~11700 K, whereas a correspond¬ 
ing MESA model at [Fe/H] = -2.20 dex, A7zahb = 0.6455 
Mq produces 10 100 K. This corresponds to a difference of 
(B — V) ~ 0.04 and a difference in ZAHB mass as recorded 
by GCB+10 of ~0.02 Mq. This difference should be in the 
correct sense, in that Pisa evolutionary tracks would record 
higher masses than the MESA models. A consistent ~0.02 Mq 
difference across all clusters would serve to increase the aver¬ 
age 77 R by «0.05 (approximately the systematic uncertainty 
we have allotted for differences among evolutionary tracks) 
but would also serve to significantly flatten both distribu¬ 
tions of 77 versus [Fe/H], It would remove the unobserved 
bright AGB stars in metal-poor clusters. However, it would 
also remove more of the bright AGB stars in metal-rich clus¬ 
ters, which are observed. 

Reducing the age of the metal-poor globular clusters 
to bring their ages in line with MF+09 and D+10 flattens 
the age distribution of globular clusters, making them much 
more closely co-eval. The addition of ~0.02 Mq to the ini¬ 
tial mass of metal-poor clusters and its subtraction from the 
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metal-rich clusters flattens the 77 distributions. The subse¬ 
quent change in 17 by ~ ±0.05 at either end is sufficient to 
remove the bright AGB stars from all but the most-metal- 
poor clusters, while making the bright AGB stars from the 
metal-rich clusters more visible. A constant age does not en¬ 
tirely solve the AGB problem but it does alleviate it, plus 
removes the unexplained gradient in both 7 ?r and 17 sc- 

We suggest that a combination of these uncertainties 
could be sufficient to remove the gradient seen in 17 , therefore 
we cannot claim any metallicity dependence exists in 17 . An 
absence of a metallicity dependence would be in keeping 
with the magneto-acoustic driving mechanism proposed for 
these winds. 

5.4-2 The spread in 17 for intermediate-metallicity clusters 

The noticeable spread in 17 for intermediate-metallicity clus¬ 
ters (-1.8 < [Fe/H] < - 1.0 dex) traces to the HB masses 
derived by GCB±10. The substantial spread in these, nearly 
0.1 Mg, is not seen in the initial masses we derive (Figure 
0. It also shows up as a significant spread in the RGB mass 
loss calculated by GCB±10 (their figure 11). 

Generally speaking, the older clusters seem to have the 
lower HB masses and higher but the underlying causes in 
this spread are unclear. Well known HB pairs (M3 and M13, 
NGC 288 and 362) lie on either side of the divide in HB mass 
at 0.635 Mg. Clusters on the lower-HB-mass, higher-?? side 
average 23 per cent higher in mass but 39 per cent sm aller in 
radius than clusters on the other side (iGnedin et al.ll2002ll . 
Higher mass and central concentration are notable charac¬ 
teristics of clusters with multiple populat ions as the y are 
more easily able to retain their ISM (iMcDonald fc Ziilstral 
l20l4l . They also have significantly larger spreads in HB 
mass (GCB+10), providing evidence of their multiple pop¬ 
ulations. The median stars in these clusters could be signif¬ 
icantly helium-enhanced, in which case our initial mass for 
the HB stars in these clusters would be over-estimated. If 
true, this would reduce the overall values of rj we derive, but 
still within the limits of the systematic errors we apply due 
to the uncertainty in helium abundance ('Appendix IC3I) . 


6 CONCLUSIONS 

In this paper, we have investigated the efficiency of mass loss 
experience by globular cluster stars and how that mass loss 
is distributed throughout the stars’ evolution. Based on the 
ZAHB masses from GCB±10, we have assumed that the 
median star experiences negligible helium enrichment. We 
have combined these isochrones with evolutionary models 
to determine the following principle conclusions: 

• Metal-rich globular cluster stars lose more mass on the 
RGB than metal-poor stars, simply due to their slower evo¬ 
lution. 

• For most giant branch stars in globular clusters, mass 
loss on the RGB dominates the total mass lost. 

• The mass-loss models of Reimers and Schroder & Cuntz 
show a relatively constant efficiency in converting (pre¬ 
sumably) magneto-acoustic energy to stellar outflow, inde¬ 
pendent of metallicity. The modelled efficiencies are 17 r = 
0.477 ± 0.070±o.° 62 and r) sc = 0.172 ± 0.0241^. The 


quoted uncertainties denote the random spread among clus¬ 
ters and the overall systematic uncertainty (including a fac¬ 
tor for helium enrichment), respectively. These values are 
towards the higher end of previous estimates in globular 
clusters, and may be somewhat higher than generally found 
in the held. 

• While a weak gradient with metallicity is possible, 
whereby metal-poor stars experience less-efficient mass loss, 
any gradient is likely to be within the systematic uncertain¬ 
ties of the underlying data. 

• Stars in globular clusters should generally lose ~0.10 
Mg in their final AGB dusty wind. Some clusters may no 
longer regularly produce any dusty stars. A more precise 
value is difficult to determine given the remaining uncertain¬ 
ties in the data. Better empirical calibration of the AGB tip 
in globular clusters is recommended. 
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APPENDIX A: INFORMATION ON 
LITERATURE STUDIES OF CLUSTER AGES 

SW02 built on methods em ployed in their earlier works 
dSalaris fc Weissl 1 19971 . 1199 81 to provide absolute ages for 
50 clusters. Four metallicity bins were chosen, with break¬ 
points on the CG97 scale at [Fe/H] = -1.75, -1.3 and -0.9 
dex. In each bin, a reference cluster was chosen (M15, M3, 
NGC 6171 and 47 Tuc) and its absolute age determined 
from the U-band magnitude difference between the red side 
of it s HB and the main-se quence turn-off (MSTO) (taken 
from iRoscnbcr g" et alJll999l h This is referred to in later lit¬ 
erature as the “vertical” method. For each bin, relative ages 
were calculated corresponding to the difference in (V — I) 
or (B — V) colour between the MSTO and the base of the 
RGB. This is referred to a s the “horizontal” met hod. Be¬ 
spoke isochrones are used dSalaris fe Weissl Il998l l. assum¬ 
ing a helium fraction of Y = 0.23 + 2>Z (for metal fraction 
Z ) and an alpha-enhancement of [a/Fe] = ±0.4 dex. This 
provides ages in good agreement with the earlier study of 
iRosenberg et alj dl999l l. which we do not include due to its 
relatively small number of clusters (35) and similar set of 
data to SW02. Ages are provided for both metallicity scales. 

DA±05 calculated ages from F439W and E5551U Hub¬ 
ble Space Telescope ( HST) and ground-base d VI photom¬ 
etry of 55 clusters dRosenberg et al.l l 2000 bl lal: IPiotto et aid 
|2002J). They follow the same approach as SW02, using five 
metallicity bins, broken at [Fe/H] = 1.8, 1.5, 1.3 and - 

1.1 dex, with templates NGC 4590, 5262, 5904, 1851 and 
6362. Minor diferences exist in measuring the ZAHB magni¬ 
tude for metal-rich clusters ([Fe/H] > 1 dex, relying on the 
lower envelope of the HB) and metal-poor clusters (relying 
on matching to a template at that metallicity). They use the 
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isochrones of IPietrinferni et al.l (l2004l i and relative ages are 
provided for both metallicity scales. 

MF+09 published ages for 64 globular clusters, on 
the basis of F606W and E814IF HST photometry 
llSaraiedini et all 120071 : lAnderson et al.l l2008l l . This is the 
same dataset and methodology which provides the space- 
based data for our horizontal branch masses, thus we invoke 
it as a primary age calibrator too. The depth of the HST ob¬ 
servations allow the authors to fit loci to the clusters’ main 
sequences up to the luminosity of the HB. By overlaying 
them in colour-magnitude space, using the MS and RGB 
base as reference points, an intrinsic difference in MSTO 
magnitude can be found. This is analagous to the “verti- 
cal” method descr ibed above. Results are compared to the 
iDotter et alJ ll2007f) isochrones fo r both metallicity scales , 
and to the evolution ary models o flPietrinferni et al.l (l2004h : 
iBertelli et al.l (Il994f) and lGirardi et ah 1 200(1 ') in the CG97 
scale. 

D+10 provide absolute ages for 61 globular clusters us¬ 
ing the same HST Advanced Camera for Surveys (ACS) 
data as MF+09. D+10 takes the ages of 55 of its clus¬ 
ters from fits to the IDotter et al.l (1200711 isochrones, which 
uses the ZW84 metallicity scale, but includes an extra six 
outer -halo glob ulars with ages calculated using the updated 
iDotter et all (I 2 OO 8 I) models (AM-1, Eridanus, NGC 2419, 
and Palomar 3, 4 and 14). Absolute ages were calibrated by 
scaling to a nominal maximum age of 13.3 Gyr. 

VBLC13 also use the HST ACS data to derive ages 
for the original 55 globular clusters. Distances are set from 
the ZAHB magnitude, while ages are derive d from isochr one 
fittin g with appropriate metal mixtures dVandenBerg et all 
I2012IT Particular care is given to the slopes of the sub-giant 
branch (SGB) and RGB in terms of the effects of chem¬ 
istry, etc., though this was found to have little effect on 
the age derived as the MSTO is largely unaffected. Their 
approach, however, can be reduced to simply adopting ver¬ 
sions of the “vertical” or “horizontal” (or both) methods 
described above. 


APPENDIX B: A COMPARISON OF STELLAR 
EVOLUTION MODELS 

B1 Existing stellar evolution models 

The choice of underlying stellar evolution model and its 
applicability to the globular clusters in question is a large 
source of systematic error in our work. In this section, we 
explore the differences between our MESA models and five 
published different sets of isochrones. We will discuss the 
uncertainties the choice of stellar evolution models provides 
in estimating the initial mass of ZAHB stars in the clusters. 

These are namely: _ 

_ The Dartmouth stellar evolution database (IDotter et al.l 

|2008| 'R. This is the original database generated for and 
calibrated against the HST ACS data mentioned above 
(MF+09, D+10). Along with the usual parameters of age 
and metallicity, it allows the user to vary the helium abun¬ 
dance between a primordial value (Y = 0.245 + 1.5Z) and 


3 http://stellar.dartmouth.edu/models/isolf_new.html 


two fixed values (Y = 0.33, 0.40) and the a-element abun¬ 
dance between [a/Fe] = -0.2 to +0.8 dex in steps of 0.2 
dex. 

The Padova database of stellar evolutionary tracks and 
isochrone From t his database, we use version 1.1 of the 
PARSEC isochrones of lBressan et al ] (ml). These isochrones 
assume a solar-scaled composition, with Y = 0.2485 +1.78Z, 
based on Zq = 0 .0152. While different options for interstel¬ 
lar and circumstellar dust extinction and carbon stars exist, 
we do not need nor implement them here. No variation of 
[a/Fe] is possible, but the PARSEC isochrones do allow the 
user to change the RGB mass-loss rate via Refiners’ 77 (see 
Section An alternative set of models is also available 
using the C Mlfl and YZV ATifl i nterfaces, from lMarigo et ahl 
d2008ll and IBertelli et all (l2008l l. respectively. A metallicity 
range from Z = 0.0001 to 0.07 is covered at Y = 0.23 to 
0.46, with an interpolation tool provided on the website at 

a user-specified value of Reimers’ 77 ._ 

_ The Pisa evol utionary library (ICastellani et al.l 120031 : 

ICariulo et al.l [2004f1. This is an older set of models a grid 
of fixed metallicities and ages, with no provided interpola¬ 
tor. Y — 0.23 is assumed for Z < 0.001; Y = 0.232 for 
Z = 0.001; Y = 0.237 and 0.27 for Z = 0.004 and Y = 0.25 
for Z = 0.008. Some variation in mixing length and convec¬ 
tive overshooting is applied to the higher-metallicity mode ls. 

The Victoria Regina models IVandenBerg et aI]|2006i 'FI. 
A variety of sub-solar evolutionary tracks calculated at 
[a/Fe] = 0.0, +0.3 and +0.6 dex, plus a set of solar-scaled 
models covering higher metallicities. An interpolation tool 
is provided, which has the ability to produce isochrones. 

A bag of stellar trac ks and isochrones (BaSTI, ver 5.0.1; 
IPietrinferni et allteOOdFf) . Used by DA+05 and MF+09. A 
set of stellar isochrones from [Fe/H] = -2.27 to +0.40 dex, 
scaled with Y = 0.245 for Z = 0 and increasing to their 
adopted solar value of Y = 0.273 for Z = 0.0198. Both 
canonical and non-canonical mixing is included for both 
solar-scaled and a-enhanced ([a/Fe] = +0.35 dex) abun¬ 
dances. An interpolation tool is provided for these models. 
Additional models are available for increased helium abun¬ 
dances (Y = 0.30, 0.35 and 0.40) and for extreme C, N, O 
and Na abundances (at solar helium abundance with limited 
models at Y = 0.28, 0.35 and 0.40). 


B2 Comparing of our MESA tracks with existing 
models 

Figure IbT 1 compares our MESA evolutionary tracks with the 
existing Dartmouth tracks. Evolutionary tracks for the PAR¬ 
SEC models are not publicly available, but the range of 
masses on the upper RGB is sufficiently small that the PAR¬ 
SEC isochrones are a reasonable substitute in this region. 

The MESA and Dartmouth isochrones follow each other 
reasonably well for most of the stellar evolution process. 


4 http://stev.oapd.inaf.it/cgi-bin/cmd 

5 http://stev.oapd.inaf.it/cgi-bin/cmd 

6 http://stev.oapd.inaf.it/YZVAR/ 

7 http://astro.df.unipi.it/SAA/PEL/Z0.html 

8 http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/ 
VictoriaReginaModels/ 

9 http://albione.oa-teramo.inaf.it/ 


© 9999 RAS, MNRAS OOO.Hfi22l 






























































A calibration of RGB mass-loss efficiency 15 


10000 


« 1000 

'( 7 ) 

o 

c 


o 
c/) 


c/) 

O 

c 


100 


10 



7000 6500 6000 5500 5000 4500 

Effective temperature (K) 


4000 


3500 3000 



Effective temperature (K) 

Figure Bl. Top panel: stellar models on the Hertzsprung-Russell diagram. The solid, red line shows a MESA track for a M = 0.9 Mg 
star with [Fe/H] = 1 dex, losing mass with tjr = 0.2. A comparative Dartmouth track with [Fe/H] = - 1 dex, [cn/Fe] = +0.4 dex, without 
mass loss, is shown as a dotted blue line. Bottom panel: a zoom into the RGB tip, showing a variety of MESA and Dartmouth tracks and 
Padova isochrones. The format of the labels is: [Fe/H], Y (helium fraction), [ct/Fe], r]\^. The effects of changing these three parameters 
can be seen by comparing the tracks. Also shown are lines of constant mass loss for the formulae of Reimers’ and SC05. Each dot on the 
MESA tracks denotes a loss of 0.01 Mg, the final dots being (from left to right) 0.75, 0.73 and 0.68 Mg. 


Minor differences exist, which are primarily due to the ex¬ 
act choices of helium con tent and detailed_elemental abun¬ 
dances of the models (see iDotter et all ()2007l ) for a detailed 
discussion of these effects). Particularly important here is 
the oxygen abundance, which accounts for the difference in 
main-sequence turn-off temperature. 

The lower panel of Figure [Bll shows the upper part of 
the RGB of selected MESA and other evolutionary models, 


revealing the sensitivity of the RGB tip to a variety of differ¬ 
ent parameters. The adoption of an a-enhanced model has 
moved our RGB tip to cooler temperatures (higher mass-loss 
rates) than the Padova models, though the proportional in¬ 
crease in helium abundance has partly compensated for this. 
The additional helium abundance compared to the Dart¬ 
mouth evolution tracks has meant that the tracks evolve 
faster and are truncated at lower luminosities. 
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The inclusion of mass loss also shifts the tracks to cooler 
temperatures, increasing the mass-loss rate, but decreasing 
the luminosity of the RGB tip. Faster mass loss also leads 
to marginally longer RGB evolution, with a difference of 
9.6 Myr between the r/ = 0.4 and 0.5 MESA models. The 
choice of r/ in the stellar models therefore has an impact on 
the calculated 77 we derive: a change in the stellar evolution 
model of Ar/n = 0.1 results in a change in the derived value 
of A? 7 r = 0.0175. It is on this basis that we chose to calculate 
models with two different values of r/n'- one with rjn = 0.4 
and one with = 0.5. 


B3 Comparing of our MESA tracks with 
observations 


Figure IB2I shows our MESA stellar evolution model applied 
to the globular cluster 47 Tucanae (NGC 104). We have cho¬ 
sen this cluster because it has one of the most homogeneous 
populations of all the populous clusters (GCB+10). The 
stellar abundances are well determined by recent publica¬ 
tions and closely follo w those we have assum e d for the MESA 
evolu t ionary tracks dGratton et al.l 120131: ICordero et al.l 


20141: Cerniauskas et al. 2014: Dobrovolskas et alJ 20141 : 
Laoenna et, al.ll2014l : iThvgesen et al.ll2014l . Johnson et al., 
submitted). Tem peratures of all stars abov e the MSTO have 
been derived in iMcDonald et al.l (|2011bl f . These tempera¬ 
tures have b een confirmed spect roscopically on the RGB 
and AGB by ICordero et all (l2014ll and Johnson et al. (sub¬ 
mitted) and should therefore have an absolute accuracy of 
<30 K between 4200 and 4700 K. This spectroscopic cali¬ 
bration ensures we have appropriately addresssed the small 
reddening correction towards the cluster and that we have 
a model with the correct chemical abundances. This allows 
us to confidently adopt parameters in MESA that fit the ob¬ 
served Hertzsprung-Russell diagram. 

The MESA tracks shown in Figure lB2l have been evolved 
for 12.119 Gyr, and the last 519 Myr of each track is shown. 
They are computed for our standard elemental mixture, 
scaled to [Fe /H] = -0.72 dex as appropriate for 47 Tuc 
(lHarrisll2010l . and references above), and adopt r/n = 0.45 
for their entire evolution. The MESA isochrone is marginally 
better at reproducing the main-sequence turnoff than the 
equivalent Dartmouth isochrone without affecting the lu¬ 
minosity of the sub-giant branch, although the accuracy of 
both is within the uncertainties imparted by cluster and 
and distance. Both isochrones reproduce the location of the 
RGB tolerably well. The MESA isochrone over-estimates the 
luminosity of the RGB bump, which is due to our si mplistic 
treatment of mixing in these stars (cf. Ca ssisi et alJ [20021. 
However, this does not have an important effect on later 
stellar evolution for the purpose of this paper. 

Further up the giant branch, at cooler temperatures, 
the Dartmouth isochrone and MESA track diverge, primarily 
due to the mass loss experienced by the MESA model. Qual¬ 
itatively, the Dartmouth model appears to better represent 
the upper RGB. However, non-LTE and dynamical effects 
become important as the stars are pulsating. Detailed treat¬ 
ment of molecular bands becomes important below 4000 K. 
Dust production also artificially lo wers the computed stella r 
temperature in the brightest stars (IMcDonald et al.l l2011a). 
As the temperature in Figure lB2l are computed with static 
atmospheres under the assumption of LTE and no dust, tem¬ 


peratures below 40 00 K may be less trust worthy. A more re¬ 
fined treatment bv lLebzelter et all (l2014l l suggests that the 
brightest, pulsating stars may have temperatures slightly 
(~100 K) lower than predicted, although there are consid¬ 
erable systematic uncertainties preventing us stating this 
with certainty. 

The post-RGB evolution is largely controlled by the 
mass-loss treatment on the RGB. Qualtitatively, our T] = 
0.45 model provides the appropriate temperature and lumi¬ 
nosity of the horizontal branch in 47 Tuc. The start of the 
AGB, at ~135 Lg is not well reproduced by the MESA track 
(80 Lg), which is largely also due to our treatment of mix¬ 
ing. The model fails to converge after three thermal pulses, 
at which point it has a mass of 0.567 Mg. A fourth thermal 
pulse could be expected, which is not expected to be enough 
to form a carbon star. Carbon stars are not observed in 47 
Tucanae. 


APPENDIX C: UNCERTAINTIES IN r) 

Calculation of the formal systematic uncertainty in ri is dif¬ 
ficult. Differences in the treatment of stellar evolution are 
important, but the main source of systematic uncertainty 
in 77 is caused by assumptions about the stellar parame¬ 
ters themselves. An exhaustive study, requiring many self- 
consistent grids of stellar evolution models, is computation¬ 
ally prohibitive. For the purposes of this study, an adequate 
approximation can be made by comparing pre-existing grids 
of stellar evolution models and test cases of MESA evolution 
models. 

Since r] is measured by following stellar evolution mod¬ 
els, comparing an initial mass with a current mass, we have 
divided our error analysis into these three categories. 


Cl Uncertainties arising from choice of evolution 
model 

Errors associated with the choice of evolution model can be 
approximated by choosing models with the same parameters 
but from different evolution codes. Figure [Cl] shows the dif¬ 
ference in initial mass of stars at the RGB tip between each 
set of isochrones and the Dartmouth models. To arrive at an 
accurate comparison, we have had to interpolate over both 
age and metallicity. We do this by two-dimensional linear in¬ 
terpolation between points in age and [Fe/H], interpolating 
past the end of existing models if necessary. Note that many 
of the vertical features in the diagrams may arise from this 
interpolation as some models are sparsely sampled in metal¬ 
licity. Where [Fe/H] is not directly given, we have calculated 
it from the metal abundance assuming Zq = 0.0152. 

The systematic error caused by the choice of stellar evo¬ 
lution model is taken to be the standard deviation of the 
RGB tip masses of all five evolutionary codes for a partic¬ 
ular age and metallicity. It is shown as the top panel in 
Figure [C2l For the age-metallicity range covered by globu¬ 
lar clusters, the error is typically ~0.01 Mg, which imparts 
a systematic error of ~0.05 to 77 . Note that exact values are 
used in the derivation of the systematic error in the main 
text. 


© 9999 RAS, MNRAS OOO.II1I221 




































A calibration of RGB mass-loss efficiency 17 




6500 6000 5500 5000 4500 4000 3500 3000 


Effective temperature (K) 

Figure B2. A Hertzsprung-Russell diagram of 47 Tucanae (grey points from iMcDonald et alifioilbl ). A 12-Gyr, Y = 0.25, [a/Fe] = 
+0.2 dex, [Fe/H] = —0.75 Dartmouth isochrone is shown as a dotted blue line. Stellar evolution tracks from MESA are shown in red. 
Tracks at several different masses are shown between 11.6 Gyr and 12.119 Gyr in order to approximate a stellar isochrone appropriate 
for 47 Tuc. From bottom (dark) to top (light), the tracks are for masses 0.83, 0.84, 0.85, 0.86, 0.867 and 0.875 Mg. 


C2 Uncertainties in initial mass 

Once the evolutionary model is set, the primary uncertain¬ 
ties arise from the stellar physics and chemistry included in 
the model. These include the helium abundance, [a/Fe] ratio 
and treatment of convective overshooting. Figure [C2l shows 
our adopted uncertainties for these effects. The quadrature 
sum of these four error components (evolution model, helium 
abundance, [a/Fe] and convective overshooting) is taken as 
the total systematic error in the age-initial-mass relation. 
This error is shown in Figure [3] in the main text (Section 
Erg. An additional error is present on rjsc that is not present 
on r /r due to the fact our model grid is calculated using 
Reimers’ mass-loss law. This is explored in more detail in 
Section IC4I 

C2.1 Helium abundance in evolution models 

The helium abundance has a relatively small effect on the 
Hertzsprung-Russell diagram, resulting in a marginally hot¬ 
ter RGB and a less-luminous RGB tip. A good comparison 
can be made between two Dartmouth isochrones at 11.5 
Gyr, [Fe/H] = 1.0 dex and [a/Fe] = +0.4 dex, for two dif¬ 
ferent helium abundances. The isochrones at Y = 0.25 and 
Y = 0.33 (AT = 0.08) show a difference in the RGB tip 
position of 40 K and 70 Lq. By far the stronger factor is 


the rate of stellar evolution. The same isochrones provide 
initial stellar masses of 0.877 and 0.750 M@ at the RGB tip. 
For a typical HB mass of 0.65 Mg, the helium-rich model 
produces an r/ that is only 44 per cent of the base model. 

This is also illustratable in the mesa models. Figure 
IC3l shows the RGB tip in the Hertzsprung-Russell diagram 
for three evolutionary tracks with marginally different he¬ 
lium abundances, each with M = 0.90 Mg. Sodium and 
magnesium abundances have been increased by 5Ay and 
oxygen and aluminium abundances decreased by the same 
amount to account for the observed Na-0 and Mg-Al anti¬ 
correlations that are thought to be associated with helium 
enrichment. The stellar parameters near the RGB are negli¬ 
gibly different for the three different cases, but a difference 
of ±0.005 in helium fraction results in a difference of /g'gg 
Gyr in stellar age and lo'ocmo ’ n ^^B tip mass. 

Since we observe clusters at a known age, rather than 
a known mass, we must convert this into a mass differential 
to understand its effect on r/. A change in the same model 
of ±0.02 Mg in initial mass produces a change in stellar age 
of 7 j q 2 Gyr. We can therefore assume that a change of AY 
= 0.005 produces a change in initial mass at fixed age of 
~ ±0.011 M 0 . 
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Figure Cl. Difference between the Dartmouth RGB tip masses 
and other isochrones, for [a/Fe] = 0 isochrones with solar-scaled 
helium abundances. The models are (top to bottom): PARSEC, 
Pisa, Victoria Regina and BaSTI. 


Figure C2. The major sources of error contributing to determi¬ 
nation of the the initial mass of RGB tip stars. Top panel: the 
error associated with choice of stellar evolution model. Second 
panel: combined random and systematic errors associated with 
the natal helium abundance. Third panel: the absolute difference 
between [cc/Fe] = +0.2 and +0.4 dex models, indicating the un¬ 
certainty due to a varying [a/Fe]. Bottom panel: from the BaSTI 
models, the error due to convective overshooting prescription. 
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Effective temperature (K) 

Figure C3. The RGB tip as a function of helium abundance 
for [Fe/H] = -0.60 dex models at 0.90 Mq with 77 # = 0.45. See 
Figure m for details of the symbols. The open squares mark a 
mass of 0.61 Mg. 


C2.2 Error introduced by the uncertain helium abundance 

GCB+10 model that most clusters have a median helium 
richness close to the primordial value (68 per cent have 
hmed ^ 0.254, compared to an assumed average of Y = 
0.247). We can therefore adopt an uncertainty in each clus¬ 
ter’s helium enhancement of AY = 0.007. As this is purely 
due to helium enrichment, not depletion, these uncertainties 
make the overall error asymmetric, reducing AM and r/. 

For our adopted error on each cluster, we use the dif¬ 
ference betwee n the Y — 0.23 and 0.26 Padova models 
dBertelli et al.ll2008l l . scaled by 0.23x as appropriate for our 
uncertainty of AY = 0.007. The resulting error is around 
0.01 Mq in most cases, rising to ~0.02 Mq for the highest- 
metallicity cluster. The corresponding values for test cases 
mentioned in the previous section are ~0.011 Mq for the 
metal-poor Dartmouth test case and ~0.015 Mq for the 
metal-rich MESA test case, i.e. very similar uncertainties are 
found among all three models. 

In addition to the individual errors, there in a substan¬ 
tial systematic error in the helium abundance attributable 
to the relative enrichment of helium with metallicity. Fol¬ 
lowing the PARSEC models, we have adopted AY/AZ = 
1.78 for our evolutionary tracks, with a primordial Y = 
0.2485. Primordial helium abundances scatter by about 
±0.001, while modern scalings of the initial solar helium 


abundance tend to vary by ±0.004 ( Asplund et al. 

200S 

Lodders et al. 20091; Serenelli & Basu 20ld: Coc et al. 

2013 

Aver et al. 20131; Planck Collaboration et al. 2013i). Mean- 


while, the same sources give Zq to be uncertain by around 
4 per cent (0.0152 ± 0.0006). The uncertainty in primordial 
helium fraction is negligible compared to other sources of 
uncertainty, but the scaling of the helium enrichment with 
Z becomes important at higher metallicities. 

We here assume that the first generation of globular 
cluster stars have experienced a helium enrichment which is 
proportionally similar to the Sun. This may not be the case, 
but the uncertainty in the relationship is sufficient to take 
into account the variations imposed by GCB±10. Based on 
the primordial and solar helium abundances, we can then 


derive an uncertainty in the helium content of these first- 
generation stars as: 

AFsystematic = 0.001 + 0.263Z. (Cl) 

Based on the calibration found in the last section, this can 
be translated to an systematic uncertainty in initial mass of: 

AMinit, systematic — 0.0022 ± 0.0090Z/Z Q . (C2) 

Following the abundances we adopted in our MESA models, 
this then becomes: 

AMinit, systematic — 0.0022 + 0.0158 x 10 [Fe/H1 . (C3) 

C2.3 Uncertainties arising from other chemical and 
physical choices 

To estimate the uncertainty introduced by assuming a fixed 
a-enhancement ([a/Fe] = ±0.3 dex), we use the difference 
between the two Dartmouth models at ±0.2 and ±0.4 dex. 
The range of ±0.1 dex is close to the standard d evia tion in 
globular clusters’ [a/Fe] of 0.11 dex llRoedieer et alJl2013h . 
The uncertainty arising from changes in [a/Fe] dominates 
the uncertainty budget for clusters around [Fe/H] = -2 dex 
and between [Fe/H] m 1 and -0.4 dex. Outside these areas, 
the choice is dominated by the choice of evolutionary model. 

Finally, we include the uncertainty due to convective 
over-shooting by taking the difference between BaSTI mod¬ 
els including two different levels of over-shooting. Again the 
effects are strongest at higher metallicities, and particularly 
become important above [Fe/H] = — 1. 

For the most-metal-poor clusters, the error in initial 
mass is set by the choice of stellar evolution model. At 
[Fe/H] ss -2 and ss - 1.2, the [a/Fe] and choice of evolu¬ 
tion model dominate. Beyond [Fe/H] = -T, the [a/Fe] and 
convective overshooting become more important and, in the 
most-metal-rich clusters, the convective overshooting and 
helium abundance dominate the errors. 

To summarise, we have adopted the following random 
uncertainties in our calculation of initial mass: 

• Age: random error following variations among literature 
derivations, as prescribed in the main text; 

• [Fe/H]: total error of ±0.1 dex (see main text), to ac¬ 
count for uncertainties in derivations for individual clusters; 

• Helium abundance (Y): fixed (/lo oi 6 Mq); 

• [a/Fe]: variable (up to ~0.01 Mq); 

and the following systematic uncertainties: 

• [Fe/H]: an error due to calibrating an absolute metal¬ 
licity scale incorporated in random error; 

• Choice of evolutionary model: variable (~0.01 Mq); 

• Helium abundance (Y): variable (0.0022 ± 0.0158 x 
10 [Fe/H1 ); 

• Convective overshooting: variable (up to ~0.02 Mq). 

C3 Uncertainties in horizontal branch mass 

Figure [C4l shows the effects of helium abundance on ZAHB 
loci and HB evolutionary models using data f rom two evo lu¬ 
tionary mod els. The PGPUC ZAHB models llValcarce et al.l 
I 2012 L l2013bh use a fixed initial mass, leading to differences 
in the core mass. The Padova synthetic HB tracks use a 
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Figure C4. Solid red lines show zero-age HB loci for stars be¬ 
tween 0.60 Mg and 0.61 Mg (increasing as shown) for differ¬ 
ent helium fractions (V). Data are from the PGPUC models for 
[Fe/H] = 1.25 dex, [ct/Fe] = +0.3 dex. Blue dotted lines show 

HB and post-HB evolutionary tracks for Z = 0.001, M = 0.60 
Mq stars from the Padova isochrones, with time increasing as 
shown. 


fixed core and total mass. For our instantaneous view of a 
globular cluster, the PGPUC models give the more accurate 
representation of where stars will land on the ZAHB, while 
the Padova models can be used as a representation of the 
directions in which those stars will evolve. 

The derivation of HB mass in GCB+10 uses purely 
colour terms (their equations (3) through (8)), which can be 
translated almost directly into effective temperatures. The 
length of the solid, red PGPUC model lines in Figure lGll rep- 
resents a 0.01 Mq change in ZAHB mass. The models are 
separated by 1.5-2.5 per cent in Y. with the lines almost 
overlapping in effective temperature. For modest increases 
in Y, up to about Y = 0.270, the lines continue to overlap. 
Thus, an increase of Y up to this value will have a <0.01 
Mq effect on the ZAHB mass derived by GCB+10. This 
increase, AY ~ 0.03, is much more than either the typical 
random variation among clusters of AY ss 0.007 found by 
GCB+10, or the projected uncertainty in helium enrichment 
with metallicity. We therefore estimate that the uncertainty 
in ZAHB mass derived by GCB+10 is <<0.01 Mq. 

An additional, second-order effect can be seen in the 
Padova models, where helium-rich stars are shown evolving 
to hotter temperatures while on the HB (Figure [C4I) . HB 
stars spend most of this period near the ZAHB locus, so the 
error in derived mass caused by changes to the HB evolution 
should be smaller than that caused by the change in ZAHB 
locus. We therefore combine both these errors to estimate 
that the ZAHB mass of the clusters as measured by GCB+10 
is ~0.01 Mq. 


C4 An additional error on r /sc 

Our evolutionary tracks are calculated based on Reimers’ 
r/. Applications to other mass-loss laws (not implemented in 
mesa) acquire an additional uncertainty. That uncertainty 
is difficult to assess properly, however it can be estimated 
using a comparison between our r /r = 0.4 and r /r = 0.5 
tracks. 



Mass at r| R = 0.4 (M 0 ) 


Figure C5. Error in mass caused by improper treatment of RGB 
mass loss. Evolutionary time runs right to left. The underlying 
evolutionary track is a model at [Fe/H] = 1, M = 0.90 Mq. The 
dashed blue line shows the mass attained using g r = 0.4, applied 
to a = 0.5 track. The solid red line shows the mass attained 
using rjsc = 0.158, applied to a tjr = 0.4 track. 


Figure [C5] shows the effect of applying the wrong value 
of r] to a track. Computing the stellar mass using g r = 0.4, 
applied to an evolutionary track computed for g r = 0.5 
results in an RGB tip mass that is 6 per cent too low. Since 
the associated change in L and R are comparatively small, 
this translates into a fractional error of 6 per cent on LR/M, 
hence also on the value of g r derived. For a derived r/ = 0.5, 
for example, the associated additional uncertainty would be 
0.03. 

Figure IC5l also shows the difference in mass found when 
applying r/sc = 0.158 to a r/B. = 0.4 model. The value of 
rjsc here is chosen to make the RGB tip masses identical. 
The difference in mass is always less than 2 per cent. Since 
the difference in L, R, T and g are comparatively small, we 
estimate that the additional fractional uncertainty associ¬ 
ated with applying the SC05 law to a track computed for 
Reimers’ law is <2 per cent. To be conservative, we adopt a 
2 per cent error in the main text. 


REFERENCES 

Anderson J., Sarajedini A., Bedin L. R., King I. R., Piotto 
G., Reid I. N., Siegel M., Majewski S. R., Paust N. E. Q., 
Aparicio A., Milone A. P., Chaboyer B., Rosenberg A., 
2008, AJ, 135, 2055 

Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, 
ARA&A, 47, 481 

Aver E., Olive K. A., Porter R. L., Skillman E. D., 2013, 
11, 17 

Bertelli G., Bressan A., Chiosi C., Fagotto F., Nasi E., 
1994, A&AS, 106, 275 

Bertelli G., Girardi L., Marigo P., Nasi E., 2008, A&A, 484, 
815 

Bladli S., Hofner S., Nowotny W., Aringer B., Eriksson K., 
2013, A&A, 553, A20 

Boyer M. L., McDonald I., van Loon J. T., Gordon K. D., 
Babler B., Block M., Bracker S., Engelbracht C., Hora 
J., Indebetouw R., Meade M., Meixner M., Misselt K., 


© 9999 RAS, MNRAS 000.1111221 





A calibration of RGB mass-loss efficiency 21 


Oliveira J. M., Sewilo M., Shiao B., Whitney B., 2009, 
ApJ, 705, 746 

Boyer M. L., Woodward C. E., van Loon J. T., Gordon 

K. D., Evans A., Gehrz R. D., Helton L. A., Polomski 

E. F., 2006, AJ, 132, 1415 

Bressan A., Marigo P., Girardi L., Salasnich B., Dal Cero 
C., Rubele S., Nanni A., 2012, MNRAS, 427, 127 
Caloi V., D’Antona F., 2007, A&A, 463, 949 
Cariulo P., Degl’Innocenti S., Castellani V., 2004, A&A, 
421, 1121 

Carney B. W., 1996, PASP, 108, 900 

Carretta E., Bragaglia A., Gratton R., D’Orazi V., Lu- 
catello S., 2009, A&A, 508, 695 
Carretta E., Gratton R. G., 1997, A&AS, 121, 95 
Cassisi S., Salaris M., Bono G., 2002, ApJ, 565, 1231 
Castellani V., Degl’Innocenti S., Marconi M., Prada Mo¬ 
roni P. G., Sestito P., 2003, A&A, 404, 645 
Catelan M., Grundahl F., Sweigart A. V., Valcarce 
A. A. R., Cortes C., 2009, ApJ, 695, L97 
Coc A., Uzan J.-P., Vangioni E., 2013, ArXiv e-prints 
Cordero M. J., Pilachowski C. A., Johnson C. I., McDonald 

1., Zijlstra A. A., Simmerer J., 2014, ApJ, 780, 94 
Cranmer S. R., Saar S. H., 2011, ApJ, 741, 54 
Dalessandro E., Salaris M., Ferraro F. R., Mucciarelli A., 

Cassisi S., 2013, MNRAS, 430, 459 
d’Antona F., Bellazzini M., Caloi V., Pecci F. F., Galleti 

5., Rood R. T., 2005, ApJ, 631, 868 

de Angeli F., Piotto G., Cassisi S., Busso G., Recio-Blanco 
A., Salaris M., Aparicio A., Rosenberg A., 2005, AJ, 130, 
116 

Dobrovolskas V., Kucinskas A., Bonifacio P., Korotin S. A., 
Steffen M., Sbordone L., Caffau E., Ludwig H.-G., Royer 

F. , Prakapavicius D., 2014, A&A, 565, A121 

Dotter A., Chaboyer B., Jevremovic D., Baron E., Ferguson 
J. W., Sarajedini A., Anderson J., 2007, AJ, 134, 376 
Dotter A., Chaboyer B., Jevremovic D., Rostov V., Baron 
E., Ferguson J. W., 2008, ApJS, 178, 89 
Dotter A., Sarajedini A., Anderson J., Aparicio A., Bedin 

L. R., Chaboyer B., Majewski S., Marin-Franch A., Milone 
A., Paust N., Piotto G., Reid I. N., Rosenberg A., Siegel 

M. , 2010, ApJ, 708, 698 

Dupree A. K., Avrett E. H., 2013, ArXiv e-prints 
Dupree A. K., Smith G. H., Strader J., 2009, AJ, 138, 1485 
Fenner Y., Campbell S., Karakas A. I., Lattanzio J. C., 
Gibson B. K., 2004, MNRAS, 353, 789 
Fusi Pecci F., Bellazzini M., 1997, in The Third Conference 
on Faint Blue Stars, Philip A. G. D., Liebert J., Saffer R., 
Hayes D. S., eds., p. 255 

Gesicki K., Zijlstra A. A., Hajduk M., Szyszka C., 2014, 
A&A, 566, A48 

Girardi L., Bressan A., Bertelli G., Chiosi C., 2000, A&AS, 
141, 371 

Gnedin O. Y., Zhao H., Pringle J. E., Fall S. M., Livio M., 
Meylan G., 2002, ApJ, 568, L23 
Gratton R. G., Carretta E., Bragaglia A., Lucatello S., 
D’Orazi V., 2010a, A&A, 517, A81 
Gratton R. G., D’Orazi V., Bragaglia A., Carretta E., Lu¬ 
catello S., 2010b, A&A, 522, A77 
Gratton R. G., Lucatello S., Sollima A., Carretta E., Bra¬ 
gaglia A., Momany Y., D’Orazi V., Cassisi S., Pietrinferni 
A., Salaris M., 2013, A&A, 549, A41 
Grevesse N., Sauval A. J., 1998, Space Science Reviews, 85, 


161 

Groenewegen M. A. T., 2014, A&A, 561, Lll 
Habgood M.-J. J., 2001, PhD thesis, The University of 
North Carolina at Chapel Hill 
Harris W. E., 1996, ApJ, 112, 1487 
Harris W. E., 2010, ArXiv e-prints 

Henyey L., Vardya M. S., Bodenheimer P., 1965, ApJ, 142, 
841 

Kalirai J. S., Saul Davis D., Richer H. B., Bergeron P., 
Catelan M., Hansen B. M. S., Rich R. M., 2009, ApJ, 
705, 408 

Karakas A., Fenner Y., Sills A., Campbell S. W., Lattanzio 
J. C., 2006, Memorie della Societa Astronomica Italiana, 
77, 858 

Lapenna E., Mucciarelli A., Lanzoni B., Rosario Ferraro 
F., Dalessandro E., origlia L., Massari D., 2014, ArXiv 
e-prints 

Lebzelter T., Nowotny W., Hinkle K. H., Hofner S., Aringer 
B., 2014, A&A, 567, A143 

Lodders K., Palme H., Gail H.-P., 2009, Landolt Bornstein, 
44 

Marigo P., Girardi L., Bressan A., Groenewegen M. A. T., 
Silva L., Granato G. L„ 2008, A&A, 482, 883 
Marin-Franch A., Aparicio A., Piotto G., Rosenberg A., 
Chaboyer B., Sarajedini A., Siegel M., Anderson J., Bedin 
L. R., Dotter A., Hempel M., King I., Majewski S., Milone 

A. P., Paust N., Reid I. N., 2009, ApJ, 694, 1498 
McDonald I., Boyer M. L., van Loon J. T., Zijlstra A. A., 

2011a, ApJ, 730, 71 

McDonald I., Boyer M. L., van Loon J. T., Zijlstra A. A., 
Hora J. L., Babler B., Block M., Gordon K., Meade M., 
Meixner M., Misselt K., Robitaille T., Sewilo M., Shiao 

B. , Whitney B., 2011b, ApJS, 193, 23 

McDonald I., Johnson C. I., Zijlstra A. A., 2011c, MNRAS, 
416, L6 

McDonald I., van Loon J. T., 2007, A&A, 476, 1261 
McDonald I., van Loon J. T., Decin L., Boyer M. L., Dupree 
A. K., Evans A., Gehrz R. D., Woodward C. E., 2009, 
MNRAS, 394, 831 

McDonald I., van Loon J. T., Dupree A. K., Boyer M. L., 
2010, MNRAS, 405, 1711 

McDonald I., van Loon J. T., Sloan G. C., Dupree A. K., 
Zijlstra A. A., Boyer M. L., Gehrz R. D., Evans A., Wood¬ 
ward C. E., Johnson C. I., 2011d, MNRAS, 417, 20 
McDonald I., Zijlstra A., 2014, ArXiv e-prints 
McDonald I., Zijlstra A. A., Boyer M. L., 2012, MNRAS, 
427, 343 

Meszaros S., Avrett E. H., Dupree A. K., 2009, AJ, 138, 
615 

Moehler S., Koester D., Zoccali M., Ferraro F. R., Heber 
U., Napiwotzki R., Renzini A., 2004, A&A, 420, 515 
Momany Y., Saviane I., Smette A., Bayo A., Girardi L., 
Marconi G., Milone A. P., Bressan A., 2012, A&A, 537, 
A2 

Monaco L., Bellazzini M., Bonifacio P., Ferraro F. R., Mar¬ 
coni G., Pancino E., Sbordone L., Zaggia S., 2005, A&A, 
441, 141 

Paxton B., Bildsten L., Dotter A., Herwig F., Lesaffre P., 
Timmes F., 2011, ApJS, 192, 3 
Paxton B., Cantiello M., Arras P., Bildsten L., Brown E. F., 
Dotter A., Mankovich C., Montgomery M. H., Stello D., 
Timmes F. X., Townsend R., 2013, ApJS, 208, 4 


© 9999 RAS, MNRAS OOO.mt22l 


22 I. McDonald, A A Zijlstra 


Pietrinferni A., Cassisi S., Salaris M., Castelli F., 2004, 
ApJ, 612, 168 

Piotto G., King I. R., Djorgovski S. G., Sosin C., Zoccali 
M., Saviane I., De Angeli F., Riello M., Recio-Blanco A., 
Rich R. M., Meylan G., Renzini A., 2002, A&A, 391, 945 
Piotto G., Villanova S., Bedin L. R., Gratton R., Cassisi 
S., Momany Y., Recio-Blanco A., Lucatello S., Anderson 
J., King I. R., Pietrinferni A., Carraro G., 2005, ApJ, 621, 
777 

Planck Collaboration, Ade P. A. R., Aghanim N., 
Armitage-Caplan C., Arnaud M., Ashdown M., Atrio- 
Barandela F., Aumont J., Baccigalupi C., Banday A. J., 
et ah, 2013, ArXiv e-prints 

Reimers D., 1975, Memoires of the Societe Royale des Sci¬ 
ences de Liege, 8, 369 

Richer H. B., Fahlman G. G., Ibata R. A., Pryor C., Bell 

R. A., Bolte M., Bond H. E., Harris W. E., Hesser J. E., 
Holland S., Ivanans N., Mandushev G., Stetson P. B., 
Wood M. A., 1997, ApJ, 484, 741 

Roediger J. C., Courteau S., Graves G., Schiavon R., 2013, 
ArXiv e-prints 

Roediger J. C., Courteau S., Graves G., Schiavon R. P., 
2014, ApJS, 210, 10 

Rosenberg A., Aparicio A., Saviane I., Piotto G., 2000a, 
A&AS, 145, 451 

Rosenberg A., Piotto G., Saviane I., Aparicio A., 2000b, 
A&AS, 144, 5 

Rosenberg A., Saviane I., Piotto G., Aparicio A., 1999, AJ, 
118, 2306 

Rutledge G. A., Hesser J. E., Stetson P. B., Mateo M., 
Simard L., Bolte M., Friel E. D., Copin Y., 1997, PASP, 
109, 883 

Salaris M., Weiss A., 1997, A&A, 327, 107 
—, 1998, A&A, 335, 943 
—, 2002, A&A, 388, 492 

Sarajedini A., Bedin L. R., Chaboyer B., Dotter A., Siegel 
M., Anderson J., Aparicio A., King I., Majewski S., 
Marfn-Franch A., Piotto G., Reid I. N., Rosenberg A., 
2007, AJ, 133, 1658 

Schroder K.-P., Cuntz M., 2005, ApJ, 630, L73 

—, 2007, A&A, 465, 593 

Serenelli A. M., Basu S., 2010, ApJ, 719, 865 

Smith G., Dupree A., Strader J., 2014, ArXiv e-prints 

Sneden C., 2004, MmSAI, 75, 267 

Thygesen A. O., Sbordone L., Andrievsky S., Korotin S., 
Yong D., Zaggia S., Ludwig H.-G., Collet R., Asplund 
M., D’Antona F., Melendez J., D’Ercole A., 2014, ArXiv 
e-prints 

Cerniauskas A., Kucinskas A., Bonifacio P., Andrievsky 

S. M., Korotin S. A., Dobrovolskas V., 2014, ArXiv e- 
prints 

Valcarce A. A. R., Catelan M., Alonso-Garcia J., Cortes 
C., De Medeiros J. R., 2013a, ArXiv e-prints 
Valcarce A. A. R., Catelan M., De Medeiros J. R., 2013b, 
A&A, 553, A62 

Valcarce A. A. R., Catelan M., Sweigart A. V., 2012, A&A, 
547, A5 

VandenBerg D. A., Bergbusch P. A., Dotter A., Ferguson 
J. W., Michaud G., Richer J., Proffitt C. R., 2012, ApJ, 
755, 15 

VandenBerg D. A., Bergbusch P. A., Dowler P. D., 2006, 
ApJS, 162, 375 


VandenBerg D. A., Brogaard K., Leaman R., Casagrande 
L., 2013, ApJ, 775, 134 

Villanova S., Geisler D., Piotto G., Gratton R. G., 2012, 
ApJ, 748, 62 

Villanova S., Piotto G., Gratton R. G., 2009, A&A, 499, 
755 

Woitke P., 2006, A&A, 460, L9 

Zinn R., West M. J., 1984, ApJS, 55, 45 


© 9999 RAS, MNRAS OOO.HU221 


